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ABSTRACT 

Context. The far-infrared (FIR) lines are important tracers of the cooling and physical conditions of the interstellar medium (ISM), and 
are rapidly becoming workhorse diagnostics for galaxies throughout the universe. There are clear indications of a different behavior 
of these lines at low metallicity that needs to be explored. 

Aims. Our goal is to explain the main differences and trends observed in the FIR line emission of dwarf galaxies compared to more 
metal-rich galaxies, and how this translates in ISM properties. 

Methods. We present Herschel PACS spectroscopic observations of the [C n] 157 /jm, [O i] 63 and 145 pm, [O hi] 88 pm, [N ii] 122 
and 205 pm, and [N hi] 57 pm fine-structure cooling lines in a sample of 48 low-metallicity star-forming galaxies of the guaranteed 
time key program Dwarf Galaxy Survey. We correlate PACS line ratios and line-to-GiR ratios with Gir, Gir/G> metallicity, and FIR 
color, and interpret the observed trends in terms of ISM conditions and phase filling factors with Cloudy radiative transfer models. 
Results. We find that the FIR lines together account for up to 3 percent of Gir and that star-forming regions dominate the overall 
emission in dwarf galaxies. Compared to metal-rich galaxies, the ratios of [O mlgg/lN ii] L 22 and [N iii] 57/[N n]i22 are high, indicative 
of hard radiation fields. In the photodissociation region (PDR), the [C n]i57/[0 1]63 ratio is slightly higher than in metal-rich galaxies, 
with a small increase with metallicity, and the [O i] i 4 s /[0 i]g3 ratio is generally lower than 0.1, demonstrating that optical depth effects 
should be small on the scales probed. The [O iii] 88 /[0i]63 ratio can be used as an indicator of the ionized gas/PDR filling factor, and 
is found to be ~4 times higher in the dwarfs than in metal-rich galaxies. The high [C ii]/Ltir, [O i]/Gir, and [O iii]/Gir ratios, which 
decrease with increasing Gir and Gir/G, are interpreted as a combination of moderate FUV fields and a low PDR covering factor. 
Harboring compact phases of a low filling factor and a large volume filling factor of diffuse gas, the ISM of low-metallicity dwarf 
galaxies has a more porous structure than that of metal-rich galaxies. 

Key words, galaxies: dwarf - infrared: ISM - ISM: H n regions, photon-dominated regions - techniques: spectroscopic - radiative 
transfer 


1. Introduction 

Low-metallicity environments, such as nearby star-forming 
dwarf galaxies, are ideal laboratories for examining conditions 
for star formation in objects as chemically unevolved as high- 
redshift galaxies. They are good candidates with which to study 
and address fundamental questions such as: How do the physical 


* Herschel is an ESA space observatory with science instruments pro¬ 
vided by European-led Principal Investigator consortia and with impor¬ 
tant participation front NASA. 


conditions of star-forming regions vary as a function of the el¬ 
emental enrichment? What controls the energy balance, heating 
and cooling, in the different interstellar medium (ISM) phases 
and shapes the morphology of dwarf galaxies? What are the roles 
of the different ISM phases in the integrated view of galaxies? 

Spectroscopy in the far-infrared (FIR) provides a unique and 
unobscured view of the ISM properties and conditions for star 
formation. The FIR cooling lines, less affected by dust attenua¬ 
tion than optical lines, are powerful probes of the star formation 
activity, linking them directly to the surrounding medium from 
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which these stars are born. The fine-structure lines [C n] 157 pm, 
[O i] 63 pm, and [O i] 145 pm are important coolants of moder¬ 
ately dense («h > 10 cm 3 ) warm (7 ~ 100 - 300 K) neutral 
ISM and are among the brightest cooling lines in star-forming 
galaxies. The [Nn] 122 and 205 pm and [Oiii] 52 and 88 pm 
lines are tracers of ionized gas. Together, they can be used as 
diagnostics of the FUV flux, gas density, temperature, and fill- 


and their general properties. We summarize the main character¬ 
istics here. The galaxies of the sample have metallicities rang¬ 
ing from ~1/4OZ 0 (IZwl8) to near solar (He2-10), and star 
formation rates ranging from ~5xl() 4 M e yr 1 (UGC4483), 
to ~25 M g yr~' (Haro 11). They have different shapes (irregu¬ 
lars, spirals, etc.), and all of them are nearby galaxies, no fur¬ 
ther than 200 Mpc. The PACS (70, 100, and 160/im) and SPIRE 


ing factor of the photodissociation regions (PDR) and ionized (250, 350, and 500 pm) photometry probe the dust emission in 

regions (e.g.. 

Tielens & Hollenbach 1985; Wolfire et al. 1990 

the FIR/submillimeter (submm). These data are essential to con- 
strain spectral energy distribution (SED) models and investigate 

Kaufman et a 

. 2006). In combination with MIR features ([S hi 


18.71 and 33.48 pm, [Nem] 15.56 pm, [Nev] 14.32/an, [Sin 
34.82 pm, H 2 , polycyclic aromatic hydrocarbons or PAHs, etc.), 
they provide a comprehensive view on conditions and excitation 
processes in star-forming regions. 

The [C 11 ] 157 pm line has been studied in a variety of ob- 


jects, such as Galactic PDRs (e.g., Bennett et al. 

1994[ |Pineda 

et al.|2013], normal spiral galaxies ( 

Stacey et al. 

1991; Malho- 

tra et ak|1997; Negishi et al. 12001; 

Brau 

ter et alJ2008), low- 

metallicity galaxies (Poglitsch et al. 

1995 

Madden et al. 1997 

Hunter et al. 2001 [Israel & MaloneyJ20L 

Cigan et al. submit- 

ted to ApJ), ultra-luminous IR galaxies (U 

JRG, Luhman et al. 

1998 

Farrah et al. 2013), and high-redshift galaxies (Maiolino 

et al. 

2009; Stacey et al. 12010; Hailey-Dunsheath et al. 2010 

Swin 

tank et al. 2012; De Breuck et al. 2014). Frc 

m those stud- 
the [Cn] line 

ies (see Genzel & Cesarsky 2000 for a review). 


emerges as the brightest cooling line of the ISM in star-forming 
galaxies and as a good tracer of the star formation activity, ac¬ 
counting for —0.1 -1 % of the FIR luminosity (e.g., Stacey et al. 
|Pierini et al.|2003] Boselli et al.||2002[ |De Looze et al. 


2010 


2011 


2014a| Sargsyan et al.|201~4 1 . A deficit in the FIR line in 


tensities is, however, found in luminous IR galaxies (e.g., Luh- 


man et al.|2003[|Gr acia-Carpio e t al.|2011| Croxall et al.|2012 


Dfaz-Santos et ak||2013 !. The origin of the [C 11 ] line has been 

debated in extragalactic studies as it can arise from molecular, 
neutral atomic, and ionized gas phases. Characterizing unam¬ 
biguously the ISM conditions thus requires the use of additional 
FIR cooling lines (such as [N 11 ] or [Oi]). In dwarf galaxies, the 
FIR lines were only detected by the KAO and ISO LWS in bright 
or Local Group dwarf galaxies (e.g., Poglitsch et al.|1995; Israel 


|et al.|1996HMadden et al.|1997||Hunter et al.|2001| i. These stud¬ 
ies found that the FIR lines are exceptionally bright compared 
to the dust emission. The brightness of the FIR lines contrasts 
with the faint CO emission, resulting in high [Cn]-to-CO(l-0) 
ratios in low-metallicity galaxies (Madden 2000] Cormier et al. 


2014| Madden et al.|in prep.|). While CO suffers from photodis¬ 


sociation, large amounts of H 2 can potentially be present in low 
metallicity PDRs owing to self-shielding. This H 2 gas not traced 
by CO is known as the CO-dark gas (e.g., Wolfire et al.|2010 i. 
Interpreting the FIR observations in terms of excitation led to the 
first evidence of such a reservoir of dark gas, not seen by CO but 
by C + ( [Poglitsch et al.|[T995 Israel et aTT|| 1 996[ 


Madden et al. 


dust and gas properties in larger samples of chemically un¬ 
evolved systems. The Herschel Guaranteed Time Key Program 
“The Dwarf Galaxy Survey” (DGS; Madden et al. 20131 is a 
230 hours project which observed 50 dwarf galaxies with the 
PACS ( jPoglitsch et al.|2010) > and SPIRE ( jClFiffin et al.|2010|> in- 


struments. Details about the program are presented in Madden 


et al. (|2013[>, which also contains a list of the dwarf galaxies 


the dust properties. The photometry of the DGS is presented in 
Remy-Ruyer et al. (2013[) and has revealed, among other results, 


1997]. This demonstrates the potential of the FIR lines as call 
brators of the dark gas at low metallicity and star formation trac¬ 
ers ( |De Looze et al.|2014a] that could be used for high-redshift 
studies. 

The advent of the Herschel Space Observatory ( |Pilbratt| 
et al. 2010 1 has provided a unique opportunity to study the 


the presence of warmer dust than in metal-rich galaxies, the pres¬ 
ence of an excess emission in the submm in several dwarf galax¬ 
ies, and lower dust-to-gas mass ratios than expected at the lowest 
metallicities ( Remy-Ruyer et al.|2014) l. The PACS (55-210 pm) 
and SPIRE ( 95-670 pm) spectrometers give access to the FIR 
fine-structure cooling lines, the CO ladder, from J — 4 - 3 to 
J - 13-12, and fainter lines of other molecules such as H 2 O, 
HCN, OH + , etc. By probing different ISM phases, they provide 
information on the filling factors and physical conditions of these 
phases, as well as on the ISM morphology. The dust and gas ob¬ 
servations together reveal a coherent multiphase picture of the 
ISM. 

In this paper, we present the observations and results of the 
FIR lines obtained with the PACS instrument onboard Herschel 
in 48 galaxies of the DGS. Section[2]describes the data reduction 
and flux extraction. The observations are interpreted empirically 
with correlation diagrams and compared to the ISO-LWS work 
of Brauher et al. ( |2Q08] ) on mostly solar metallicity objects in 
Seek [3] The data are further analyzed with H 11 region/PDR mod¬ 
els in Sect.[4]to contrast the physical conditions characteristic of 
the ISM of dwarf galaxies to those of more metal-rich galaxies. 
We discuss our results in terms of ISM properties and filling fac¬ 
tors of the ionized and neutral phases in Sect. [5] More specific 
modeling of the PDR conditions in the dwarf galaxies will be 
performed in a subsequent paper that we refer to as Paper II. 

2. Data 

2.1. Herschel PACS spectroscopy 

2.1.1. Observing details 

The PACS spectrometer observed 48 galaxies as part of the DGS 
and SHINING programs, in a total time of 152.7 h. Only the two, 
faint galaxies Tol 0618-402 and UGCA 20, which are part of the 
DGS sample (50 galaxies), were dropped from the PACS spec¬ 
troscopy program because of time constraints. The PACS array 
consists or 5x5 spatial pixels (or spaxels) covering a total field of 
view of 47" x47". The beam size is ~9" and 12" and the spec¬ 
tral resolution 90kms and 240 km s 1 at 60 pm and 160 pm, 
respectively (Poglitsch et al .|2010[ ). We obtained pointed obser¬ 
vations for the most compact objects, small mappings for more 
extended galaxies, and partial mappings of specific star-forming 
regions for the most extended galaxies {dis >6' typically), such 
as the Magellanic Clouds. We define an extended source sam¬ 
ple for the galaxies partially mapped or for which the size of 
the maps do not match for all observed lines. This sample com¬ 
prises 4 Local Group galaxies, IC 10, NGC 6822, LMC (8 re¬ 
gions targeted), and SMC, and the nearby galaxy NGC 4449. The 
other 43 galaxies of the DGS are classified as compact (see also 
Sect. [2X7]. 

The dwarf galaxies were observed between November 2009 
and August 2012. 46 sources were done in chop-nod mode and 
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6 sources in unchopped mode. The number of targeted spectral 
line varies from just [Cn] 157 pm, in the faintest cases, up to 

7 lines, [Cn] 157 pm, [Om] 88 pm, [Oi] 63 pm, [Oi] 145 pun, 
[Niii] 57 pm, [Nil] 122yum, and [Nn] 205 pm, in the brightest 
galaxies. We provide details on the observations, such as OB- 
SIDs, map sizes, lines observed, etc. in Appendix [A] 


2.1.2. Properties of the FIR lines 


General characteristics of the FIR fine-structure cooling lines ob¬ 
served with PACS are given in Table [T] 


[Cu] 157 pan: 

The [C n] line is one of the most important coolants of the ISM, 
as carbon is the fourth most abundant element. [C n] is excited 
by collisions with e~, hydrogen atoms, or molecules. The ion¬ 
ization potential of C° being 11.26eV, below that of hydrogen, 
C + can be found outside of Hii regions, in the neutral phase. It 
requires only 91.3 K to be excited hence it can cool any warm 
neutral phase. Thus C + can originate from diffuse ionized gas 
as well as diffuse neutral gas or the surface layers of PDRs (up 
to Ay ~5 mag). The relative contribution of each medium to its 
observed intensity is a function of density and ionization degree 
(e.g., Kaufman et al.[2006) . 


[Nn] 122 pan and [Nn] 205 pm: 

[N n] is only found in the ionized gas. The critical densities with 
e - are 310 cnT 3 and 50 cm -3 for [N n] 122 pm and [N n] 205 pm, 
respectively. Being in the same ionization stage, their ratio is a 
good electron density tracer of the low-density diffuse ionized 
gas. Because [C ii] can be found in the low-density ionized gas, 
part of its emission can correlate with the [N ii] emission. Hence 
the [N ii] lines can, in principle, be used to disentangle the frac¬ 
tion of [C ii] from the ionized gas to that from the PDR alone 
( |Oberst et al.|2006) . 


[Nm] 57pm: 

The [N m] 57 pm emission only arises in ionized gas and is usu¬ 
ally associated with Hu regions. Being at two different ioniza¬ 
tion stages, the ratio of [N iii]57/[N ii] 122 is a measure of the ef¬ 
fective temperature of the ionizing stars (|Rubin et al.|1994|. 


[0 1 ] 63pm and [Oi] 145pm: 

0° has an ionization potential of 13.62 eV, just above that of hy¬ 
drogen. The first two fine-structure transitions require excitation 
energies of 228 and 325 K above the ground state, correspond¬ 
ing to the [Oi] 63 pm and 145 pm lines. [Oi] is only found in 
neutral gas and usually arises from warm, dense regions. Along 
with [C 11 ], [O i] 63 pm is one of the brightest PDR cooling lines 
(e.g., [Bernard-Salas et al.||2012] >. [Oi] can also exist deeper in 
the PDR than [C 11 ] since the formation of CO occurs at larger 
Ay (~ 10 mag) than the transition of C + into C° (~3mag). The 
ratio of the two [0 1 ] lines is an indicator of the gas temperature 
for temperatures in the vicinity of 300 K. It is a density tracer for 
high temperatures and high densities. The [Oi] emission can be 
affected by optical depth effects (self-absorption, optical thick¬ 
ness), particularly for the lower level transition at 63 pm (Liseau 


et al. 2006|l. In the optically thin limit, the [0 1 ] 63 pm, just above 


the ground state, is brighter than the 145 pm line. 


[O hi ] 88 pm: 

The [O in] 88 pm line is only found in the ionized gas, and be¬ 
cause it requires energetic photons (35 eV), it is commonly ac¬ 
cepted that it comes from Hii regions rather than inter-clump 
diffuse media where radiation fields are expected to be softer. 
The [Oiii] transition at 52 pm (' Pi - 'Pi ) is at the edge of the 
PACS wavelength coverage and was not observed in our survey. 


Table 1. Characteristics of the PACS FIR fine-structure cooling lines. 


Species 

A 

[pm] 

Transition 

IP 

[eV] 

A E/P a) 

[K] 

^crit 

[cnT 3 ] 

[C ii] 

157.7 

~P 3/2 - ~ P 1/2 

11.3 

91 

50 (i) , 2.8 x 10 3 

[Nii] 

121.9 

3 P 2 - 3 P, 

14.5 

188 

310 

[Nn] 

205.2 

2 ° 

1 

k o 

O 

14.5 

70 

48 

[Niii] 

57.3 

3 Py2 - 3 P 1/2 

29.6 

251 

3.0 x 10 3 

[Oi] 

63.2 

3 Pi - 3 p 2 

- 

228 

4.7 x 10 5 

[Oi] 

145.5 

3 Po - 3 Pi 

- 

327 

9.5 x 10 4 

[Oiii] 

88.4 

3 P| - 3 Po 

35.1 

163 

510 


Values taken from |Madden et al.| ( |2013| l. The IP for [On] is 
13.62 eV. (a) Excitation temperature AE/k required to populate the 
transition level from the ground state. ( b ) Critical density for collisions 
with electrons. 


2.1.3. Data reduction 

The data were downloaded from the Herschel Science Archive 
and reduced with the Herschel Interactive Processing Environ¬ 
ment (HIPE, |Qtt||2010| > user release 12 and the PACS calibra¬ 
tion tree v65. We reprocessed the observations from the raw 
data (level 0) to spectral data cubes (level 2) using the standard 
pipeline scripts available within HIPE. We opted for the tele¬ 
scope background normalization calibration method for all chop- 
nod observations, except for the [N ii ]205 observations which 
were processed with the de-leaked RSRF provided in HIPE0 In 
the first part of the pipeline (level 0 to level 0.5), the astrometry 
and information from the instrument are read. Masks for bad and 
noisy pixels, and for when the chopper and grating are moving 
are also applied. From level 0.5 to level 1, glitches are flagged, 
and the calibration is applied. In the chop-nod case, signal from 
the chopper position is subtracted. From level 1 to level 2, a 
spectral flat-held correction is applied, and outlier data points 
are masked by a sigma-clipping. The flat-field scaling is both 
multiplicative and additive for chop-nod observations, and ad¬ 
ditive only for unchopped observations. In the unchopped case, 
a transient correction and background subtraction are also exe¬ 
cuted. For the transient correction, we applied a multi-resolution 
method as described in |Lebouteiller et ak| ( |2012| ). For the back¬ 
ground subtraction, we used a procedure that smoothes and av¬ 
erages the OFF-spectra taken before and after an ON-spectrum 
in time (when several OFF-frames are available), and then sub¬ 
tracts it from the ON-frame. 


2.1.4. Spectral maps 

The data are saved as spectral cubes before they are rebinned and 
projected onto a final spatial grid in the pipeline script. Once 
the spectral cubes (one cube per raster position) are exported 
from HIPE, we used the software PACSman version 3.57 for the 
analysis (Lebouteiller et al .j2012| ). 

The first step is to fit the spectral lines and measure the line 
flux in all the spatial positions. A polynomial curve of order 1 
or 2 is first fit to the continuum, and then a simultaneous fit 
to the continuum and line is performed with the IDL procedure 
mpfitfun ( |Markwardt|2009 ). The line profile is assumed to be 
Gaussian, although skewed profiles can be seen if the observed 
object is a point-source not centered on the array. Only UM 448 


more information available in the PACS Observer's Manual (v.2.5.1) 
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shows clear asymmetric and wider profiles in several spaxels, 
but since we have made mini-maps for this galaxy, the better 
centered spaxels are well-enough fitted by a Gaussian profile, 
and we can reliably recover the total flux of this galaxy. By de¬ 
fault, the continuum is defined as the part of the spectrum within 
2-6 FWHM of the line center. The width of the Gaussian is 
taken as the instrumental width, plus a broadening when neces¬ 
sary. The PACS line profiles are sometimes found larger than the 
instrumental profile, necessitating broadened line fits. We report 
on the broadening parameters of the FIR lines in Appendix [E] 

To estimate the uncertainty on the fit parameters, we use a 
Monte Carlo approach. We randomly perturb the data within 
their error bars, fit the line, iterate the process 100 times, and take 
as final fit parameters the median of the resulting values, and as 
error on those parameters the standard deviation. An additional 
15% absolute flux calibration uncertainty is to be considered 1 . 

For the irregular galaxy IC 10, the [Cn] spectra are con¬ 
taminated by emission from the Milky Way. This emission 
appears at +320 km s 1 relative to IC 10, with a peak of 
~1.8Jy, and is fairly uniform across our [Cn] map. The sep¬ 
aration between the two emission peaks is large enough to 
clearly identify the respective contributions. Thus to extract 
the flux only from IC 10, we fit simultaneously two Gaus- 
sians. We also note that OFF-spectra are not available for sev¬ 
eral observations in LMC-30Doradus (OBSIDs 1342231280, 
1342231282, 1342231283, 1342231285), LMC-N11B (OB¬ 
SIDs 1342188940, 1342188941), and NGC4449 (OBSIDs 
1342197813, 1342197814), resulting in variations across wave¬ 
length of a few Jy typically in the baselines of the ON-spectra. 
For NGC 4449 and LMC-N11 B, we model those variations with 
a sinusoidal component of low frequency to improve the baseline 
fit. For LMC-30 Dor, no sinusoid is required as those variations 
are not significant compared to the line strength and thus do not 
affect our line fits. 

Once flux values are calculated for each spatial position, the 
data are projected on a new grid of sub-pixel resolution of 3.1" 
(roughly 3 times oversampled) with a drizzling scheme. When 
more than one spaxel falls on the same sub-pixel, a fraction of the 
flux from the corresponding spaxel of each raster is considered, 
and the uncertainties are combined quadratically. The final flux 
maps, with spectra and line fits are shown for each galaxy in 
Appendix [B] 

2.1.5. Flux extraction 



o o o o z z z 


Fig. 1. For each PACS spectral line, the number of sources observed and 
detected (hashed) in the Herschel Dwarf Galaxy Survey. 

of V9 to account for the number of sub-pixels per independent 
beam. We consider a spectral line detected if the S/N is >3, and 
upper limits are given as the 3cr level, from the method used for 
the detected lines or from the central method otherwise. Table [4] 
lists line fluxes as well as the method used for the flux extraction 
(aperture size when applicable). To summarize, these values are 
the most reliable measurements of the total line fluxes for the 
DGS galaxies. For the extended sample, we report peak surface 
brightnesses instead of total fluxes. 

2.1.6. Summary of detections 

The [C it] 157 line was targeted in all objects, along with the 
[O hi] 88 and [O i]g 3 lines for almost all objects (~85%). For the 
brightest galaxies (30%), the [O i] 145 and [N 11 ] 122 lines were also 
observed, and in a very few cases (7%) the [N 11)205 and [N 111)57 
lines as well. Figure |T| shows a histogram of how many galax¬ 
ies were observed and detected (>3cr level) in each spectral line. 
The [C 11 ] 157, [0 111 ] 8g, and [0 1)53 lines are the brightest, detected 
in more than 90% of the cases, and on average are at least 10 
times brighter than the other lines. The [N 11 ] 205 pm line has 
only been attempted a few times but not detected as it is located 
at the edge of the spectral band, suffering from leakage. 


We have explored several methods to extract accurate line fluxes. 
For compact objects, we consider 4 methods to measure the total 
line flux, as it depends on how well-centered the source is on the 
PACS array and on the spatial extent of the source: 

- central: we apply a point-source correction only to the central 
(or brightest) spaxel. This method works only if the source is 
well centered on the spaxel. 

- psf: we fit a point spread function (PSF) to the array. 

- 3x3: we combine the inner 3x3 spaxels before fitting the line 
and apply a 3x3 point-source correction. 

-5x5: we combine the 5x5 spaxels of the PACS array before 
fitting the line and apply a 5x5 point-source correction. 

For final flux, we select the method that gives the highest re¬ 
sponse with a good signal-to-noise ratio. 

In the case of extended or elongated objects we simply perform 
an aperture extraction with sky subtraction when possible or con¬ 
sider the total flux in the area mapped. The errors are added 
quadratically from the error maps and multiplied by a factor 
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2.1.7. Additional sources of uncertainty 

In addition to the measurement errors and uncertainties from the 
PACS calibration, there are other sources of uncertainty to con¬ 
sider when interpreting the data. We discuss two main factors. 

Inhomogeneous line mapping: 

Although we have measured the total flux within each line map, 
the spatial coverages of those maps, in a given galaxy, may not 
match exactly. To estimate whether the total line fluxes can be di¬ 
rectly compared with each other, we consider the 100 pm PACS 
photometry maps as a reference for the distribution of the FIR 
emission. For each spectral line, we measure the fraction of the 
100 pm emission seen inside the map coverage. We find that this 
fraction does not vary by more than 15% for all galaxies, except 
for NGC 1140 for which we consider a correction factor of 20% 
because the [C 11 ] 157 map covers a larger area than the other lines. 
Hence we can confidently use the total fluxes reported in Table[4] 
to analyze line ratios. 
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Incomplete mapping with PACS: 

For 11 galaxies of the compact sample, the [C n] 157 line is de¬ 
tected at a >5<x level at the edges of the PACS map. The [C 11 ] 157 
emission is extended but not fully mapped, hence we may be 
missing signal when reporting total fluxes. We coarsely esti¬ 
mate the flux loss by measuring the dust emission from the 
PACS 100 pm maps inside and outside of the area covered by the 
[C 11 ] 157 maps. This way, we find a missing flux fraction of 15% 
for Mrk 153, Mrk 1089, Mrk930, NGC 1569, and of 55%, 35%, 
35%, 20%, and 60% for NGC 2366, NGC 4214, NGC 4861, 
NGC 625, and UM 311, respectively, and a fraction that we con¬ 
sider negligible (lower than 15%) for the other galaxies. Table[4] 
reports the uncorrected fluxes, and the correction factor to apply 
if one wants to compare the line fluxes to integrated quantities. 

2.1.8. Comparison with ISO 

Three galaxies of the DGS (Haro 11, NGC 1569, and NGC 625) 
were detected by ISO and reported as compact in Brauher et al. 
(2008). The LWS and PACS fluxes in Haro 11 agree within 20- 
40%, the [Oi ]63 and [O iii]xx fluxes in NGC 1569 agree within 
10%, while the [C 11)157 fluxes are found higher with PACS by 
a factor of 2 and 3 in NGC 625 and NGC 1569, respectively. 
In those galaxies (especially NGC 1569), the PACS maps show 
clear extended structure, beyond the LWS beam. Since they were 
considered as point-sources in the LWS flux extraction and the 
LWS beam profile is more peaked at 160 pm than at the shorter 
wavelengths, this may explain why the [C 11)157 fluxes recovered 
by PACS are higher. Given that mini maps were performed with 
PACS for those three galaxies, we are confident that we have 
better recovered their total fluxes. More generally, no systematic 
differences between the LWS and PACS fluxes have been estab¬ 
lished by the PACS instrument team. Therefore in Sect. [3] we 
compare our PACS data with the LWS data of other galaxies as 
they are reported in the literature (without correction factors). 


2.2. Ancillary data 

2.2.1. Spitzer IRS spectroscopy 

For the modeling carried out in Sect. [4] we compiled mid-IR line 
fluxes for the compact sources of the DGS, measured both with 
the low-resolution modules [R ~ 60 - 127) and with the high- 
resolution modules (R as 600) of the IRS instrument ( |Houck| 
et al. ]2004[ ). The line fits were performed assuming a Gaussian 
profile and constraining the line width to be the instrumental 
broadening since there is no evidence of lines being resolved. 
The IRS line fluxes are provided in Appendix |C| We provide in 
the following some details on the spectral extraction. 

The low-resolution data were retrieved from the CASSIS 
spectral database v 6 (Lebouteiller et al. 20l I [ ). For point-sources, 
we used the optimal extraction spectra, which provide the best 
S/N. For compact sources, we used the tapered column extrac¬ 
tion, in which the flux is integrated within a spatial window 
whose width scales with the source extent and with wavelength. 
While scaling the width to account for all the emission along 
the long (cross-dispersion) axis allows a better flux calibration 
for extended sources, it does not correct for the emission out¬ 
side the slit along the short (dispersion) axis. For this reason, we 
have applied a custom flux calibration correction in which we 
have assumed that sources can be modeled as 2D Gaussians with 
cropped wings (representing the source function) convolved with 
the instrument PSF. The fraction of light falling outside the slit 
is then calculated for all wavelengths and for any source extent. 


High resolution spectra were extracted with two methods 
used by the CASSIS high-resolution pipeline (Lebouteiller et al. 
in prep.| l. For point-like sources we used an optimal extrac¬ 
tion similar to the low-resolution algorithm while for extended 
sources we used a full aperture extraction (integrating the flux 
falling inside the aperture). For the full-aperture extraction, we 
decided to use a point-source calibration because in all cases the 
source is dominated by emission that is not uniform across the 
slit and that is instead characterized either by a PSF profile or 
by a broad Gaussian-like profile. The full aperture extraction 
ensures that all the light entering the aperture is accounted for. 
While we should in principle also correct for the emission out¬ 
side the aperture, like for the low-resolution data, the custom flux 
calibration for the high-resolution spectra of partially-extended 
sources is not available yet. We chose therefore a correction fac¬ 
tor similar to what is used for point-sources in the lack of a better 
estimate. 

Fluxes measured in the low- or high-resolution spectra can 
differ. While the high-resolution allows in general a more reli¬ 
able determination of the line fluxes, the low-resolution module 
slits are better adapted for extended sources. More light enters 
the low-resolution long slits (SL, LL) than the high-resolution 
apertures (SH, LH), so when using a point-like source calibra¬ 
tion, the high-resolution line fluxes will be smaller than the low- 
resolution line fluxes. The more extended the source is, the larger 
this effect. Since the size ratio SH/SL (0.24) is larger than LH/LL 
(0.14), the effect is also larger for the long-wavelength modules. 
In Sect. [4] we use line ratios from high-resolution data when 
available and from low-resolution data otherwise. Line ratios 
of [Siv] 10.5//m/[S hi] 18.7 pm, [Ne in] 15.6/im/[Neii] 12.8 pm, 
and [Sin] 18.7jum/[S hi] 33.5 pm agree within ~30% between 
the low-resolution and high-resolution data. 


2.2.2. Photometry and reference sample data 

In this paper, the dwarf galaxy sample is primarily compared 
to the ISO LWS extragalactic sample of FIR fine-structure lines 
from |Brauher et al.| ( |2008[ ) (hereafter B08), which is one of 
the most complete in the literature. The B08 sample contains 
a larger variety of galaxies than the DGS, with normal star¬ 
forming galaxies, starbursts, and active galactic nuclei (AGN). 
We only use their data for unresolved sources, and with clear de¬ 
tections (>3cr). The DGS sample extends the analysis of Brauher 


et al. ( |2008| > to lower metallicities and less luminous galaxies, 
as the B08 sample is mostly composed of bright objects. Ad¬ 
ditionally, we contrast our objects with the ULIRG samples of 
the He rsche l SHINING (jGracia-Carpio etal.|201l| and HERUS 
( |Farrah et al.|2013| > surveys. 

We study the behavior of the FIR cooling lines of the dwarf 
galaxies as a function of four important parameters: metallic- 
ity, F60/F100, Ltir, and Ltir/Lb. The broadband ratio of IRAS 
flux densities at 60 pm and 1 00 yam, denoted F60/F100, probes 
the peak of the SED (which is often at wavelengths shorter than 
ID pm in our dwarf galaxies; Remy-Ruyer et al. |2013 ), and thus 
the dust temperature. Ltir is the total infrared luminosity, mea¬ 
sured from 3 to 1 100 pm. It reflects a large part of a galaxy’s 
energy budget and is an indirect probe of the star formation ac¬ 
tivity. Ljir is commonly used as the proxy for the gas heating 
(e.g., |Rubin et al.||2009] ) since Ltir includes all contributions 
from the dust emission, and in particular the small grains and 
PAHs which are the main agents of the gas photoelectric heat¬ 
ing. The ratio of the [C 11 ] and [0 1 ] lines to Ltir is a measure of 
the photoelectric efficiency, which we discuss in more detail in 
Sect. 1 5.2 1 Finally, Lb is the blue band luminosity, corrected for 
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Galactic extinction using the Schlegel et al.| ( fl998| values pro¬ 
vided by the NASA/IPAC Infrared Science Archive but not for 
internal dust attenuation, and converted from the magnitudes re¬ 
ported in |Madden et al.| ( |2013| >. Ltir/L r relates the stellar light 
processed by dust to that escaping, i.e. the mean extinction to¬ 
wards the star-forming region. The mean extinction is sensitive 
to the compactness of the cold surrounding medium and thus to 
the covering factor of PDRs. We note that we opt for the B-band 
luminosity rather than FUV in order to have a sufficient number 
of galaxies with reliable data. 

For the normal galaxies (B08, SHINING, and HERUS sam¬ 
ples), the IRAS photometry is used to compute F60/F100 and to 
calculate L TIR using the |Dale & Helouj ( |2002| ) formula. For the 
compact sample of the DGS, we consider SED integrated values 
from Remy-Ruyer et al. ( in prep.|i fo r Ltir and observed IRAS 
fluxes from Engelbracht et al. l 2008| > for F60/F100 (when avail¬ 
able). For the galaxies undetected or not observed by IRAS, we 
create synthetic IRAS fluxes based on the Herschel photome¬ 
try and SED modeling of those galaxies, by convolving the SED 
models with the IRAS response curves. For the extended sample, 
we create Ltir and synthetic IRAS maps based on the model of 
Galliano et al. ([2011[) that we convolve to the largest beam (12" 


corresponding to the PACS 160 pm beam size). For the compar¬ 
ison to the spectral lines, we integrate those quantities over the 
same area mapped for each PACS line. 

Although less sensitive, IRAS being an all-sky survey, we 
prefer to use the IRAS photometry for our analysis with FIR 
colors rather than the Herschel photometry to ease the compar¬ 
ison with the literature. We list the IRAS and Herschel /PACS 
band ratios and Ltir in Appendix [D] We find median ratios of 
IRAS60/PACS70 = 0.86 and IRAS 100/PACS100= 1.05. Be¬ 
cause the method used to derive Ltir for the B08 galaxies and 
for the dwarf galaxies is not the same (IRAS formula versus SED 
modeling), we assess its influence by comparing the results for 
those galaxies where we can apply both methods (i.e. the dwarf 
galaxy sample). We find that the difference in the resulting Ltir 
is small, with a median difference of 15%. 


3. Observed ISM properties 

In this section, we present and analyze the correlations between 
several FIR tracers to identify observational trends in the sam¬ 
ple, and link them to general parameters (e.g., metallicity, dust 
temperature). We focus our analysis on the integrated emis¬ 
sion of galaxies. To quantify possible correlations, we com¬ 
pute the Spearman’s rank correlation with the IDL procedure 
r_ correlate only for the compact sample, when it contains 
more than 8 data pairs. We identify correlations at a significance 
level of 5%. The interpretation of identified trends is supported 
by radiative transfer modeling of the photoionized and PDR gas 
in Sect. [4] 


3.1. Luminosity of the FIR lines 

The [O hi] 88 pm line is found to be very bright in the DGS 
galaxies and is the brightest line in most cases. This is true, lo¬ 
cally in extended sources where the peaks of emission are ~ 10 
times brighter than the [Cn] 157 pm line (e.g., LMC-30Dor), 
and even on galaxy-wide scales. This is surprising, at first, as it 
requires 35 eV photons to create a 0 ++ ion while the PDR lines 
are more easily excited. Bright [Oiii] is likely related to the star¬ 
bursting nature of the dwarfs, which host stars with higher effec¬ 
tive temperatures than normal galaxies, as observed by Hunter 
et al.]( |200l[ ). The [Cn]i 57 line is the second brightest FIR line 
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and [Oi ]63 is the third brightest line. [Oi ]63 usually dominates 
the cooling in resolved PDRs (e.g., the Orion Bar,|Bernard-Salas 


et al.|| 2012 i, and we find that it is often brighter than [Cn] 


157 


on peaks of emission within extended sources (e.g., SMC-N 66 ), 
but fainter than [C 11(157 on galaxy-wide scales. 


We compare the FIR line luminosities with Ltir in Fig- [2] 
The emission from the FIR lines is strongly correlated with Ltir 
(p — 0.97 in the dwarfs), a result expected since all those quan¬ 
tities scale with the size, i.e. luminosity, of the galaxy. Those 
correlations being tight further indicates that the bulk of the FIR 
emission is related to star formation. Calibrating the FIR lines 
as star formation rate tracers is done in [De Looze et al.| ( |2014al >. 
The dwarfs extend the relation observed, at high luminosities 
with the ULIRGs from HERUS and SHINING and moderate lu¬ 
minosities with the B08 galaxies, to much lower luminosities, 
thanks to the sensitivity of Herschel. The dispersion is smaller in 
the star-forming galaxies than in the ULIRGs, which are known 
to present a line deficit, especially visible in the HERUS ionic 
line fluxes ( |Gracia-Carpio et al.|2011||Farrah et al.|2013| l. What 
is striking in the figure is the offset (~ 1 dex) between the dwarf 
galaxies and the ULIRGs from HERUS and SHINING regarding 
their [O m] xh emission. The B08 galaxies fill the gap between the 
two samples. Moreover, we notice an offset of the DGS sample 
relative to the B08 sample, with fainter [N 11 ] 122 emission in the 
dwarfs. 


3.2. Effects of metallicity 

We correlate the FIR measurements with metallicity in Fig. [3] 
As expected from the metallicity-luminosity relationship, the lu¬ 
minosity of the FIR lines increases with metallicity (p ~ 0.65), 
with large scatter (dispersion ~0.8dex). For reference, we mea¬ 
sure the following slope: L(line) oc (12 + log O/H) 2 . No clear 
trends are observed between the FIR line ratios or line-to-LTiR 
ratios and metallicity, except a slight decrease of [Oi] 63 /Ltir 
(p - -0.5) and [O iii] 88 /[C 11 ] 157 (p = -0.4) with increasing 
metallicity. No trend is detected if we include the regions of 
the LMC (orange stars). The [O i] 63 /[C 11 ] 157 ratio is also found 
higher on average at lower metallicity (c.f. IZw 18, SBS0335- 
052), but the scatter is important. To first order, variations in 
metal abundances (caused for example by delayed injection of 
carbon in the gas phase compared to oxygen) are not respon¬ 
sible for the trends with [O i ] f ,3 as they would affect the other 
lines, and in particular the [0 1 ] 145 /IG 11 ] 157 ratio, similarly. Op¬ 
tical depths effects on the [O i]g 3 line (see Sect. |2.1.2| could re¬ 
duce the amount of observed [O i ](,3 emission and account for 
the trends observed with [O i]® but not with [Oi]i 45 . However, 
they should also translate into a decrease of [O i] 63 /[C 11 ] 157 with 
Ltir, which we do not observe, and we do not have enough ob¬ 
servations of the [ 01 ] 145 line at the lower metallicities to quan¬ 
tify those effects. Therefore, the weak trends that we observe 
with metallicity are likely due to a change of cloud conditions, 
i.e. increased [Oi ]63 emission in denser/warmer clumps at lower 
metallicities (where cosmic rays or soft X-rays may become im¬ 
portant heating sources; e.g., |Pequignot|2008j ). This also agrees 
with the fact that the [O i] 63 /[C 11 ] 157 and [Oi] 63 /Ltir ratios are 
higher in the extended galaxies, particularly in the Magellanic 
Clouds observations that provide a zoomed in view on the star¬ 
forming regions and for which [O i ]63 is the main coolant of 
the dense PDR. We further test those hypotheses with radiative 
transfer models in Paper II. Overall, this indicates that metallic¬ 
ity is not the only regulator of the FIR line emission. 















































Cormier et al.: PACS spectroscopy of the Herschel Dwarf Galaxy Survey 




10 '° - 
10 9 - 

" 10 8 - 
_l 

00 10 7 - 
00 

g; 10 6 - 

_l 

io 5 - 
10 " - 
io'° - 
io 9 - 
" io 8 - 

_l 

P io 7 - 

in 

i io 6 - 

_l 

io 5 - 
io" - 


p = 0.97 


Sj,. 1 ' # 

■ 

'/'i 

»*• * 


,JT 


1/50 1/30 1/20 1/9 1/3 1 Z 0 


L(TIR) [L e ] 


10 6 


10 8 10 u 
L(TIR) [L 0 ] 


10 l4 


Fig. 2. PACS line fluxes versus Ltir. The DGS galaxies are represented by filled circles (compact sample) color-coded by metallicity. The red 
asterisks are the galaxies from Brauher et al. (20081 and the blue asterisks are the low-metallicity sources originally published by Hunter et al. 
([2001 . The orange rectangles are the HERUS data from |Farrah et al.| ( |2013| > and the red rectangles are the SHINING data from Gracia-Carpio 
|et al.| 201 1ft . Spearman correlation coefficients applied to the DGS compact sample are indicated in the top left corner (when found above the 
chosen significance level and with more than 8 data points). 


3.3. PACS line ratios 

We compare the PACS line ratios with FIR color, L TIR , and 
L-yin/Ln for the DGS and B08 samples (Fig. [4]). Median, min¬ 
imum, and maximum values of the PACS line ratios are sum¬ 
marized in Table [2] We indicate values for the complete sample 
(compact and extended objects), compact sample alone, and for 
the B08 dataset. Overall, the data clearly indicate large real vari¬ 
ations in the line ratios as well as systematic offsets in the ion¬ 
ized gas tracers. Our dwarf galaxies show weak trends within the 
sample, but behave differently than the more metal-rich galaxies 
as a group (not with respect to metallicity/abundances). 


3.3.1. [Oiii]88 / [N 11 ] 122 

[0 111 ] 88 and [N u| 122 are both tracers of the ionized gas, with 
similar critical densities (~ 300-500 cm -3 ), but different ioniza¬ 
tion potentials (35.5 and 14.5 eV, respectively). It is remarkable 
how the dwarf galaxies occupy a completely different parame¬ 
ter space in [O nijss/tN n]i 22 than the more metal-rich galaxies 
of the B08 sample. The [O nijss/tN 11 ] 122 ratio is observed to be 
high in the dwarfs, with a median value of 86, which is ~30 times 
higher than in B08. Our preferred interpretation for this high ra¬ 
tio is the presence of harder radiation fields in the dwarfs. This is 
supported by the fact that for those dwarf galaxies where we have 
observed [N m], the line is also enhanced compared to [N 11 ]. Sec¬ 
ondly, the harder radiation fields naturally lead to warmer FIR 


color (high F60/F100), as is observed. [O mJss/tN 11)122 is also 
sensitive to variations in the N/O abundance ratio, which could 
be low at low metallicities. However, we find no trend of this 
ratio with metallicity (although [N 11 ] 122 was not observed in the 
lowest metallicity galaxies). 

We do not find any statistically significant trend within the 
ensemble of dwarf galaxies between [O m]8s/[N 11 ] 122 and the pa¬ 
rameters considered. With higher F60/F100 and lower L TIR on 
average, the dwarfs do not invalidate the weak trends observed 
in Brauher et al. |( ]2Q08[ >, of [O m]/[N 11 ] increasing with FIR color 
(Fig. |4jv) and decreasing with L TIR and L TIR /L B . The fact that the 
[O iii]88 and [N 11 ] 122 trends with FIR color are weak may indicate 
a combined effect of harder radiation fields and lower densities, 
with a contribution of low-density (high and low-excitation) ion¬ 
ized gas to these lines, in addition to the compact H 11 regions. 


3.3.2. [Niii ]57 / [N 11 ] 122 

The [N in] 57 pm line was observed in only four DGS sources 
(Haro 11, Haro 3, He 2-10, and LMC-N1 IB), and three of them 
have both [Nm] and [N 11 ] 122 detections. The [Niii] 57 /[Nii]i 22 
ratios are high, with an average value of 8.1, which is 4 times 
higher than in B08 (Table [2ji. This ratio is not sensitive to 
variations in elemental abundances and an indicator of the ef¬ 
fective stellar temperature (Rubin et al. 11994| >, as it requires 
14.5 eV to create N + , and 29.6eV to create N ++ . Thus higher 
[Niii] 57 /[Nii]i 22 ratios show the presence of harder radiation 
fields and are expected in dwarfs compared to normal galaxies. 
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Fig. 3. PACS line fluxes, ratios, and line-to-Zvn R ratios as a function of metallicity. The DGS galaxies are represented by filled circles (compact 
sample) and stars (extended sample), color-coded by metallicity. Spearman correlation coefficients applied to the DGS compact sample are 
indicated in the top left comer (when found above the chosen significance level and with more than 8 data points). The red, horizontal lines show 
the mean values of the |Brauher et al.| l [2008| l sample, at metallicity around solar. 


The [N hi]57/[0 hi] 88 ratio is on average 6 times lower in the 
DGS than in B08. This ratio is sensitive to ionization degree, 
N/O abundance ratio, and density in the range 10 2 - 10 4 cm 5 
(see Table [T}. Considering the high [Niii] 57 /[Nii]i 22 values, the 
data again suggest a high filling factor of low-density high- 
excitation material in dwarf galaxies. 


density and radiation field strength ( Kaufman et al.|2006| , as the 
critical density and excitation energy of [ 01 163 are higher than 
for [Cn], and therefore higher ratios are expected in compact 
star-forming regions (e.g., SBS 0335-052), which have higher 
F60/F100 and Ttir/Tb ( |Malhotra et al.|200l] l. 


3.3.3. [Oi ] 63 / [C n ] 157 

[Oi ]63 and [Cn]i 57 both trace the neutral ISM, and the ratio of 
[O i] 63 /[C n] 157 in the DGS galaxies falls in the range observed 
by B08, with a median value of 0.6. [O i] 63 /[C ii] 157 is not corre¬ 
lated with Ltir, but we see an increase of [O i] 63 /[C 11 ] 157 with 
FIR color (f> = 0.68), and a weak correlation with Ltir/Tb 
( p = 0.36) for the dwarfs (Fig. |4j?, [4]:/ ). Theoretical PDR studies 
have demonstrated that the [O i] 63 /[C 11 ] 157 ratio increases with 


3.3.4. [Oi]i45 /[Oi] 63 

The [0 1 ] line ratio is mainly sensitive to the density and tempera¬ 


ture of the neutral gas, and to the optical depth of the cloud (Tie- 
|lens & Hollenbach||1985[ ). The median ratio of [0 1 ]i4s/[0 1]63 
is 0.07 in the dwarfs. The range of values is narrow and a lit¬ 
tle higher than the median value found in the more metal-rich 
galaxies. 

We do not find statistically significant trends between the 
[0 1 ] i45/[0 i ]63 ratio and F60/F100, L T ir. or L T ir/L b (Fig. [4]?), 
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Table 2. PACS line and line-to-L TIR median ratios of the DGS galaxies. 



Full DGS sample 

Compact sample® 1 ' 

B08 sample 

PACS line ratio 

[O I] 63 /[C Il]l 57 . 

romiss/rc ni, S7 . 

0.59“' (0.15 dex) 
2.00™ (0.34 dex) 
0 -074pMi (0.06 dex) 
2.96™ (0.24 dex) 
86.3*y (0.23 dex) 
8.06“® (0.13 dex) 
0.025““ (0.23 dex) 
0.17““ (0.52 dex) 

< 1 

0.65 (0.20 dex) 
1.79 (0.34 dex) 
0.074 (0.05 dex) 
2.66 (0.22 dex) 
61.5 (0.22 dex) 
8.06 (-) 

0.022 (0.27 dex) 
0.17 (0.37 dex) 

0.72 (0.25 dex) 
0.54 (0.45 dex) 
0.063 (0.25 dex) 
0.74 (0.25 dex) 
3.27 (0.37 dex) 
1.91 (0.57 dex) 
0.12 (0.21 dex) 
1.08 (0.44 dex) 

[O 1 ] 14s/[0 i]63 . 

[O m] 88 /[0 1] 63 . 

rOmWrNnl,,,. 

[N iii] 57 /[N Il] 12 2 . 

[Nii] 122 /[Cii ] 157 . 

[N iri] 57 /[0 m] 8 8 . 

[Nll]205/[Nll]i22 . 

line-to-Ly i R ratio [in %] 

[C Il]l57/LriR . 

[Oi]63/TtIR . 

[Om] 88 /Z/nR . 

[Ol]l45/LxiR . 

[Nii]i22/LtiR . 

[N III]57/LtIR . 

([Cn ] 157 +[Oi] 63 )/Gt,r . 

([C Ii] 157 +[Ol] 6 3 +[OlIl]g 8 )/L TIR 

0.25° ® (0.24 dex) 
0.17“® (0.14 dex) 
0.50' “ (0.29 dex) 
0.011™ 3 (0.19 dex) 
0.008Ippo*8 (°- 21 dex ) 
0.042pl 4 (0.31 dex) 
0.47' “ (0.15 dex) 
1.03™ (0.18 dex) 

0.25 (0.22 dex) 
0.17 (0.14 dex) 
0.47 (0.26 dex) 
0.010 (0.12 dex) 
0.0060 (0.25 dex) 
0.038 (0.07 dex) 
0.43 (0.17 dex) 
0.97 (0.20 dex) 

0.13 (0.33 dex) 
0.10 (0.21 dex) 
0.080 (0.30 dex) 
0.0065 (0.18 dex) 
0.016 (0.19 dex) 
0.043 (0.17 dex) 
0.26 (0.26 dex) 
0.35 (0.26 dex) 


Median ratios (considering >3 cr detections) and their dispersions in parenthesis, measured as 1.5 times the median absolute deviation in logarithmic 
space. The upper and lower script values correspond to the maximum and minimum values, respectively. 

(a) DGS sample excluding the LMC and SMC regions, NGC 4449 , NGC 6822 , and IC 10 . 


contrary to Hunter et al. ( |200l| who noted a slight increase with 
L tir /L b (on a smaller dataset). The lack of correlation with FIR 
color is expected as the [O i] line ratio stays constant for a con¬ 
stant Go/«h (provided Go < 10 4 ). Ratios higher than 0.1 are 
only found in two galaxies (0.17 in VII Zw 403, at a 3<x level, and 
0.11 in LMC-N159), and are usually explained by noticeable op¬ 
tical depth effects on the [O i] 63 p m line (|Malhotra et al.|2001[ 
|Hunter et al.|2001||Rollig et al.|2006}|Kaufman et al.|2006| l. On 

integrated scales, optical depths are not high in the dwarfs (Ay 
<3 mag; |Abel et al.|2007| ). 


3.3.5. [Om] gg / [O i ] 63 


As observed in earlier studies of dwarf irregular galaxies by 
Hunter et al. (|2001|l, the [O irr] K8 /[0 r]f ,3 ratios are high in the 


dwarf galaxies, with median value 3. While Brauher et al. ( |2008| > 
noted a slight decrease of [O iii]88/[0i ]63 with Ltir/Lb, we find 
no correlation with FIR band ratio, Ljir. or Ltir/Lb within the 
ensemble of dwarf galaxies. The DGS galaxies as a group oc¬ 
cupy a different parameter space than the B08 galaxies (Fig. [^0- 
The [O iii]gg/[0 1](,3 ratio, which is insensitive to variations of 
elemental abundances, may be used as an indicator of the filling 
factor of ionized gas, traced by [O hi], relative to the filling factor 
of PDRs, traced by [O i]. In the extended sources, [Oiii] is de¬ 
tected over large spatial scales (e.g., LMC-N11B, LMC-30Dor, 
NGC 4449; Lebouteiller et al. 201 2| |Chevance et al. in prep.| 
|Karczewski et al. in prep.| ) indicating the presence of low density 
channels where the hard UV photons can travel far (>20pc) to 
excite the gas. The dwarf galaxies are probably mostly filled, by 
volume, with ionized gas, with densities likely below ~510cm 3 
(»crit °f [O hi] 88), and low filling factor of PDRs. 


3.3.6. [N n] 122 / [C II] 157 

The [N n] i 22 /[C n] 157 ratio is on average 0.02, 5 times lower 
than in the B08 galaxies. This is mostly a consequence of the 


high ionization fractions which favor the presence of N ++ rather 
than N + . The [N n] 122 /TC 11 ] 157 ratio seems uncorrelated with FIR 
color, Ltir, or Ltir/Tb in the dwarfs (Fig. [4] ). :Brauher et al. 
( |2008[ > also found an absence of correlation for all morphologi¬ 
cal types. 

This ratio can be used as an indicator of the [C 11 ] excitation 
mechanism, where high [N ii]/[C 11] ratios can mean that a non- 
negligible fraction of the [Cn] emission originates in the ion¬ 
ized gas and low ratios indicate that [C 11 ] originates in the PDR 
gas. [Oberst et al. ( |2006 1 show that ionized gas alone yields at 
least [N 11 ] 122 /tC 11 ] 157 >0.1. The conditions for [C 11 ] emission to 
originate in the ionized gas are low densities (« cr it,e— 50 cm~ 3 ) 
and low ionization, as carbon would otherwise be doubly ion¬ 
ized (I.P. 24.5 eV), i.e. the [N m] emitting gas does not pro¬ 
duce notable [Cn] emission. While Brauher et al. ( 2008] ) sug¬ 
gest that a significant fraction of the [Cn] emission in metal- 
rich galaxies originates in H 11 regions where [N 11 ] comes from, 
our lower [N 11 ]r 22 /[C 11 ] 157 values favor a non-Hu region origin 
for the [C 11 ] line, from the neutral and/or very low-density ion¬ 
ized gas. From individual analyses, we have found that: (1) for 
Haro 11, where [N n] r 22 /[C 11 ] 157 = 0.05 is the highest ratio in the 
DGS galaxies, the low-density diffuse ionized gas contributes up 
to ~40% to the [C 11 ] emission; (2) in the LMC regions, where 
[N 11 ] i 22 /[C 11 ] 157 - 0.03, [CII] is associated with the ionized gas 
at <15% (|Lebouteiller et al.||2012| |Chevance et al.||in prep.|l; 
(3) in NGC 4449, [N 11 ] 122 /[C n] 157 a 0.02 and the fraction of 
[Cn] arising from the ionized gas is small (Karc zewski et al.frn 
|prep.| l. Therefore we estimate that, typically, less than 15% of 
the [C 11 ] emission arises in the ionized gas. Although the ion¬ 
ized gas fills a large volume in the dwarf galaxies, its excitation 
state is too high for C + to be present. Lower N/O abundance at 
lower metallicity (e.g., [Liang e t ak||2006| may change this es¬ 
timate, but the fraction should remain low as we observe low 
[N 11 ] i 22 /[C ii] 157 ratios even at half-solar metallicity. [N 11 ] 122 is 
not detected for metallicities 12+log (O/H) < 7.8 and, where de¬ 
tected, [N 11 ] i 22 /[C ii] 157 does not show any trend with metallicity. 
More precise quantification of this fraction in the DGS galaxies 
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Fig. 4. PACS line ratios against F60/F100 and Ltir/Lb. The DGS galaxies are represented by filled circles (compact sample) and stars (extended 
sample), color-coded by metallicity. See also the caption of Fig. [2] 


requires multiphase modeling which is beyond the scope of this cate that star-forming regions clearly dominate the overall emis- 
paper. sion of the galaxy. 


3.4. FIR fine-structure line-to-L T i R ratios 
3.4.1. General remarks 

In Fig. [5] we show ratios of PACS line-to-L T1R versus F60/F100, 
CtifQ and Ltir/Lb to investigate variations in the cooling by the 
gas and dust. 

The line-to-L TIR ratios are found high in the dwarf galax¬ 
ies, at least twice as high as in the B08 galaxies for [C ii]/Ltir 
and [Oiii]/Ltir (see Table [2]). Among the highest points are the 
extended galaxies, in which we have zoomed on the most ac¬ 
tive star-forming regions. Excitation by the young massive stars 
results in bright FIR lines, while the dust can still emit signifi¬ 
cantly further from the active regions. Line-to-L TIR ratios would 
probably be lower for the LMC in its entirety. In more quiescent 
regions of the extended galaxies, line-to-LjiR ratios are lower in 
general. This is especially the case for [Ora] and [Oi], and less 
for [C n] because of the contribution from diffuse emission to the 
total [Cn] flux. |Rubin et al.| ( |20091 ) estimate [Cii]/Ltir -0.45% 
for the entire LMC, which is close to our average extended sam¬ 
ple value. High ratios for the compact objects may therefore indi- 


2 Even though we plot line over Ltir versus Ltir, these diagrams are 
still useful as the x-axis is an indicator of galaxy size while the y-axis 
probes the heating and cooling physics of the ISM. 
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3.4.2. [Oiii]88 / Ltir 


[Oiii]88/Ttir values span an order of magnitude, from 0.1 to 1.7, 
with average value of 0.50%. The dwarf galaxies are at the high 
end of the B08 sample. The dwarfs also extend the trend ob¬ 
served in B08 of [Oiii]88/Ttir decreasing with increasing Ltir 
and Ltir/Lb, which is weaker for the dwarfs (p - -0.4) than 
for the B08 galaxies (Fig. |5J. This suggests enhanced [Om] 
emitting volume when the escape fractions are high. However, 
there are multiple possible interpretations of those correlations. 
[O iii]/L T i R could also be reduced by a significant fraction of the 
Ltir emission unrelated to star formation. However similar be¬ 
havior is observed in [Om]/[Oi] (see Fig. [4]'). Because both of 
these lines are directly related to star formation, we rule out this 
interpretation. It is also possible that Ltir/Lb does not directly 
measure the covering factor of the star-forming regions if there 
is a significant contribution from field stars to the B-band lu¬ 
minosity. We have tested this possibility by making similar dia¬ 
grams for galaxies with available K-band (field star dominated) 
and FUV (young stars) data. We find that [Oiii]/Ltir is actu¬ 
ally more strongly correlated with Ltir/Lfuv than with Ltir/Lb 
and uncorrelated with Ltir/Lk, which favors an interpretation in 
terms of UV escape fraction. Thus the correlations are best ex¬ 
plained by the combination of high stellar effective temperature 
on galaxy-wide scales, resulting in high [O hi] emission, and en- 
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Fig. 5. PACS line-to-inR ratios as a function of F60/F100, Ctir. an d Ctir/Lb- Same legend as in the caption of Fig. [5] 


hanced ISM porosity with lower Ltir and Ltir/^b. which in turn 
allows UV photons to excite 0 ++ on larger spatial scales. On the 
other hand, [O iiiJss/^tir is uncorrelated with FIR band ratios 
in the dwarfs, while B08 find an increase of [O iii]88/^tir with 


F60/F100. They attribute this increase in warmer galaxies to a 
higher density of Hii regions. As discussed in Sect. 3.3 [Oiii] 


may also originate from low-density ionized gas, hence flatten¬ 
ing the relation in dwarfs. 
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3.4.3. [Cn] 157 / Ttir and [Oi]g 3 / Ltir 


The ratio of [C ii]/L T i R spans more than an order of magnitude 
in the DGS galaxies, from 0.04 to 0.7%, with a median value 
of 0.25%. The ratio of [Oi] 63 /Ltir does not vary significantly 
among the dwarfs, with a median value of 0.17%. Both ratios 
have median values slightly higher in the dwarfs than in B08. 
Together, the median value of ([Cn]i 57 +[Oi] 63 )/Ltir is 0.47%. 
Since [Cn] and [Oi ]63 are the main coolants of the PDR, where 
the gas heating is dominated by photoelectric effect on dust 
grains, which cool through Ltir, the ([Cn] +[0 i]63)/Ltir ratio 
is often used as a proxy for the photoelectric efficiency (e.g., 
Tielens & Hollenbach 1985 1 . The larger spread in the [C ii]/Ltir 
values compared to [Oi] 63 /Ttir (Fig. [5]» may hint towards differ¬ 
ent origins of the [C n] line (ionized/neutral gas) or different con¬ 
ditions of the [C n]-emitting phases (as [Cn] saturates for den¬ 
sities > 10 4 cirr 3 ). Since the spread in ([O i ]63 +[Cii])/Ttir is 
smaller compared to that seen in [C ii]/L T ir (0.15 dex compared 
to 0.26 dex), a mixture of dense and diffuse phases overall com¬ 
pensate for the differences in the [O i ]63 and [C n] emission. 

The only, weak, trend identified in the dwarfs with FIR color 
is a decrease of [Cii]/L T ir (p = -0.47); this has been ob¬ 
served in the B08 galaxies, although the trend is less pronounced. 
This can be interpreted as a compactness effect with reduced 
emission in denser and/or warmer environments (e.g., |Malho-| 
tra et a l.pOOl) . Moreover, [Oi] 63 /Ltir is weakly anti-correlated 
with Ltir. The dwarfs seem to extend the correlation observed in 
B08 of [C h] 157/FriR and [Oi^/L-tir decreasing with increasing 
Ltir/Tb (p — -0.62). The decrease in the [OiIm/Ltir values is 
small compared to the decrease in the [C ii]/Ltir values. As a re¬ 
sult, ([C 11 ] 157 +[0 i] 63 )/Ljir a lso decreases with increasing L TIR 
(p =; -0.52) and Ljir/Tb (p — -0.81). We discuss those high 
line-to-LxiR and correlations further in Sect. |5.2| 


3.4.4. ([C 11 ] +[Oi] 63 +[0 in]) /Ltir 

The ratio of ([C 11 ] +[Oi ]63 +[Oiii])/Ltir spans values from 0.4 
to 2.1, with a median of 1.03. This is about 3 times higher than 
in the B08 sample. The sum of [C 11 ] +[0 1 ]® +[0 in] may be re¬ 
garded as the budget of the gas cooling (not a proxy for the pho¬ 
toelectric efficiency), which is clearly enhanced in the dwarfs, 
and Ltir as the total dust cooling. 

Following the previous remarks, this ratio is uncorrelated 
with FIR color, as [O m] dominates the FIR line cooling and no 
relation was found for [Oiii]/Ltir, but it clearly anti-correlates 
with L tir and L TIR /L B (p - -0.6, -0.7). The DGS galaxies ex¬ 
tend the trends observed in the B08 sample. 


3.5. Summary of correlation analysis 

In summary, the dwarfs present high [Om]/[Nn], [Cii]/[Nii], 
[Oiii]/[Oi], [Cii]/L tir , [Oi]/L tir , and [Oiii]/L tir ratios, which 
set them apart as a group from metal-rich galaxies. Within the 
ensemble of dwarfs, we do not find evident correlations between 
their FIR emission and the considered parameters (metallicity, 
F60/F100, L tir , Lt IR /L b ), which indicates that their observed 
properties are not controlled by a single parameter. 

The interpretation of the observed trends is not straightfor¬ 
ward since several factors act and compete in the emission of the 
cooling lines and Ltir. These factors include cloud physical con¬ 
ditions, mixing of ISM phases, geometry, dust properties, and 
photoelectric efficiency. They all play a role to some extent since 
the observed ratios are values integrated over full-size galaxies. 
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4. Trend analysis with radiative transfer models 


In this section, we characterize the range of physical conditions 
describing the ISM, more specifically the Hii region and the 
PDR, of the DGS galaxies. We aim to interpret the observed 
group behavior in Sect.[3]with radiative transfer modeling, and in 
particular, to identify how the parameter space of average phys¬ 
ical conditions change compared to metal-rich galaxies. 

We note that there is no study on large samples of galaxies 
that self-consistently models the Hii region and the PDR, be¬ 
sides the work on IR luminous galaxies by Abel et al. ( 2009[) and 
Gracia-Carpio et al.| ( |20lT| ) and on individual objects in Kaufman 
et al. ( 2006| ). We aim for a general comparison of the metal-rich 
and metal-poor datasets, and leave the detailed ISM modeling of 
the dwarf galaxies for Paper II. 


4.1. Line ratios and typical conditions in metal-rich galaxies 

For the PDR properties of “normal” metal-rich galaxies, we only 
use entries from |Brauher et al.| ( |2008| ) that have L T i R between 
10 9 and 10 12 L o (~ 130 galaxies) to exclude extreme sources. 
More than half of that sample, comprising normal, starburst, 
and AGN galaxies was analyzed by Malhotra et al. ( |200 1 [ 1 and 

conditions 


in the PDR using the 
n h — IO 2 - 4 5 cm~ 3 anc 

Kaufman et al. 

(1999) models and found 

G 0 - 10 / 4:1 . In the comparison to PDR 

models, Malhotra et a 

. (2001) and Negishi et al. (2001) correct 


for a non-PDR contribution to the [C 11 ] emission which is on the 
order of 50%. In the dwarfs, we have estimated in Sect. |3.3.6| 
that this fraction is likely <15%. We apply a correction of 10% 
for non-PDR [C 11 ] emission to our sample. 

For the H 11 region properties of metal-rich galaxies, we refer 
to the studies on MIR ionic lines of the Spitzer Infrared Nearby 
Galaxies Survey (SINGS; |Dale et al.|2009)> and of the luminous 
IR galaxies from GOALS i lnami et al.|2013) . We consider neon 
and sulphur measurements for the galaxies in the luminosity 
range L T ir - 10 9 - 10 12 L Q (~190 sources). The [S m] is/[S 111)33 
ratio is traditionally used as a tracer of the electron density (e.g., 
|Dudik et al.|2007| >, and the [Ne 111 ] is/[Ne n] 12 and [S iv] io/[S 111 ] ^ 
are diagnostics of the hardness of the radiation field ( |Giveon| 
et al.|2002 Verma et al.|2003 Groves et al.|2 008). A typical so¬ 


lution for the H 11 region of those metal-rich galaxies has an elec¬ 
tron density n e - 350 cm -3 and ionization parameter log U <* — 3 
( Dale et al.|2009 Inami et al.|2013| l. For comparable constraints 
on the H 11 region properties in the dwarf galaxies, we use the 
IRS fluxes extracted for the compact DGS galaxies in Sect. |2.2| 
Table [3]lists the median line ratios for the metal-rich and metal- 
poor galaxies. 

We refer to the H 11 region and PDR conditions for normal 
galaxies as the metal-rich fiducial model (Table [3j. 


4.2. Model parameters 

We model the H 11 region and PDR simultaneously with the spec¬ 
tral synthesis code Cloudy v. 13.03 (Ferland et a l.|2013[ >. The ge¬ 
ometry is ID spherical, consisting of the central source of en¬ 
ergy surrounded by a cloud of material. The source of radiation 
is chosen as a young starburst that we simulate with Starburst99 
. Following |Dale et~akj ( |2009] >, we adopt a 
OMyr continuous star formation scenario (with a Rroupa IMF 
and Geneva stellar tracks), and fix its luminosity to 10 9 L G . 


(Leitherer et al.|201()) 
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Table 3. Observed and modeled properties of the dwarfs and normal 
galaxies. 


IRS line ratio 

Metal-rich sample 

Compact DGS sample 

[SmMSm],, 

1.43 (0.14 dex) 

1.53 (0.12 dex) 

[Si V ] 10 /[Sm],8 

0.18 (0.38 dex) 

1.30 (0.42 dex) 

[Ne m] 15 /[Neii] 12 

0.17 (0.35 dex) 

4.90 (0.48 dex) 

Metallicity 

2 Zism 

0.25 Zism 

Reference ISM abundances (Z ISM ): 


O/H = 3.2 x 10 -4 

, C/H = 1.4 x 10 -4 

, N/H = 8 x 10 -5 

Ne/H = 1.2 x 10- 

\S/H = 3.2 x lO - 

5 

H ii region fiducial model 

density [cm 3 J 

10 25 

10 2 ° 

log U 

-3.0 

-2.5 

Rs.eff [pc] 

16 

46 

PDR fiducial model 

density [cm 3 ] 

10 35 

10 40 

log G 0 

3.3 

2.7 


Median IRS line ratios with their dispersions in parenthesis, measured 
as 1.5 times the median absolute deviation in logarithmic space ( |Dale| 
|et al.|2009| and |Inami et al,|2013| for the metal-rich galaxies; this work 
for the low-metallicity, DGS galaxies). ISM abundances are taken from 
|Cowie & Songaiia| (1986| ); |Savage & Sembach| ]1996| >; |Meyer et al. 


(|1998 . The model parameters of the metal-rich sample are from 
et al.| 2009 [ and Inami et aI7|(|2013|> for the H n region, and from 


hotra et 


al. (200 


effective radfus of the 


Dale 


[MIL 


and|Negishi et al.|(2001) for the PDR. Rs.eS is the 
:Hii 


region. 


We vary the following physical conditions: the hydrogen 
density («h; from 10 10 to 10 50 cm -3 in increments of 0.5 dex) 
and the inner radius (rj n ; from 10 19 ' 7 to 10 22 0 cm, or 15 pc to 
3,200pc, in increments of 0.3 dex), which is the distance be¬ 
tween the starburst and the surface of the modeled ISM shell. 
For each model, we compute U. Go, and the IRAS F60/F100 
band ratio. The FIR color is strongly correlated with the ioniza¬ 
tion parameter U which is varied through the density and radius. 

It is defined as U = —^——, where 2(H) is the number of 
4nrf n n H c 


ionizing photons and c is the speed of light. Go is the intensity 
of the FUV radiation field taken at the PDR front, in units of 
1.6 x 10~ 3 ‘ 


erg cm 


(Habing 


1968). 


The gas-phase abundances are representative of the ISM 
(averages of observations in our galaxy), with Zjsm = 12 + 
log (O/H) = 8.5, and given in Table [3] The grain composition 
and size distribution are assumed typical of the galactic ISM, 
with Ry = 3.1 ( jWeingartner & Draine||2001| . The dust-to-gas 
mass ratio is 6.3 x 1 0" 3 and the PAH abun dance (number o f 
grains per hydrogen) is PAH/H = 1.5 x 10 -7 (Abel et al. 2008 I, 
corresponding in mass to M(PAH)/M(dust) = 0.02. The mean 


metallicity measured in the DGS sample with the Pilyugin & 
|Thuan| ( |2505l ) calibration is 12 + log (O/H) = 7.9. In the SI NGS 


galaxies, it is 12 + log (O/H) = 8.4 with the Pilyugin & Thuan 
(|2005|) cali bration and 12 +log (O/H) - 9.0 with the |Kobulnicky 


|& Kewley| ( |2004| l calibration (Mousta kas et al.|2010), with l arger 
discrepancies in the methods for massive galaxies ( Kewley & El-| 
lison|2008 1. Therefore we take for representative metallicity of 


the metal-rich galaxies 12 + log (O/H) = 8.7. We run two grids 
of models, one at metallicity 0.25 Zism for the dwarf galaxies, 
and one at metallicity 2 Zism for the metal-rich galaxies. For the 
model grids, the metallicity of the stellar tracks as well as the 
metal, grain, and PAH abundances in Cloudy are scaled to the 
metallicity of the model (0.25 or 2). We also include a turbu¬ 


lent velocity of 1.5 km s , as done in 


Kaufman et al. 


(jl999). 


All models are stopped at an Ay of 10 mag (corresponding to 
a column density of /V(H ) = 10 23 and 10 22 cm 2 for the low 
and high-metallicity grids, respectively), to ensure that they go 
deep enough in the PDR. We note that the main PDR tracers 
are mostly formed at lower Ay and changing the column density 
around those values does not change their emission significantly. 
The grids are generated assuming constant density. 

4.3. Model results 

Figure[6]shows the model results for the high and low-metallicity 
grids. To extract information on the ISM properties, we con¬ 
centrate on the following ratios: [S m]i8/[S m] 33 , [S iv] io/[S in] jg, 
[Ne hi] i 5 /[Ne n] i2 , [O m] 8S /[N n] i22 , [O iiiJss/Ttir, for the H n re¬ 
gion; [Oi]63/[Cn]i57, [Cii]i 57 /L T ir, [Oi]63/Ttir for the PDR. 
The IRAS F60/F100 band ratio and the Go of the models are 
also plotted. To ease the comparison with the literature, we plot 
quantities in the H ii region as a function of U and in the PDR 
as a function of Gq. First of all, we notice that all observed ra¬ 
tios (except [O iii]/Ltir) are covered by the model grids, hence 
we can identify a parameter space that is favored for the dwarfs. 
Fiducial model parameters for the high and low-metallicity cases 
are given in Table [3] 


4.3.1. Conditions in the H ii region: high-metallicity case 


In the high-metallicity case (left panel of Fig. |5J, the fiducial 
model for the H ii region agrees partially with the observed IRS 
line ratios of metal-rich galaxies. The [S m] 33 /[S in 1 1 8 ratio is rel¬ 
atively independent of U and the observed values agree best with 
densities between 100 and 1,000 cm -3 . Both [Nem]i 5 /[Nen]i 2 
and [S iv]io/[S m]i8 ratios are quite sensitive to U, the age of 
the burst, and less to density. The [Ne m] 15 /[Ne n]| 2 ratio is 
slightly over-predicted by the fiducial model. The median and 
dispersion values of this ratio reported in Table [3] are driven by 
the GOALS galaxies, while a few SINGS galaxies clearly have 
higher [Nem]i 5 /[Ne ii]i 2 ratios ( [Dale et al.||2009] i. |Inami et al. 
( 2013| l have investigated the influence of metallicity and star- 
formation history (continuous versus bursty) on the predicted 
[Ne in] i 5 /[Ne n] i 2 ratio. They demonstrate that their observed 
[Ne in] i 5 /[Ne n] i 2 values are best match by models with young 
instantaneous bursts rather than continuous star formation mod¬ 
els which systematically over-predict this ratio. The largest dis¬ 
crepancy between observations and model predictions is seen 
for the [S iv]!o/[S m]is ratio, which is under-predicted by one 
order of magnitude. Part of this discrepancy is due to an ob¬ 
servational bias because the detection rate of the [S iv]io line is 
lower (about 4 times less than the neon lines). For the galax¬ 
ies where [S iv]io/[S hi] i 8 is detected, the [Ne hi] is/[Ne ii] i 2 ratios 
are higher than the reported median value. The [Ne hi] is/[Ne n] i 2 
and [S iv]iq/[S m]i8 ratios are well correlated, and extrapolating 
this correlation to all galaxies where [Nem]i 5 /[Nen]i 2 is de¬ 
tected would lower the observed median ratio of [S iv] io/[S iii]is 
by a factor of ~3. Despite this correction, it seems difficult to re¬ 
produce the [Nem]i 5 /[Neii]i 2 and [S iv]io/[S in] is ratios simul¬ 
taneously with our high-metallicity grid. 

Concerning the PACS lines, both [O iii]/Ltir and 
[O m]/[N ii] i 22 are well matched by the fiducial model. This 
suggest that [Oiii]s8/[N ii]i 22 may be used as an indicator of the 
hardness of the radiation field as an alternative to the IRS ratios. 
This model also predicts that about 10% of the [Cn] emission is 
originating in the ionized gas, while Malhotra et al. (2001]) and 
Negishi et al.| (|2001|) estimated a ~50% contribution using the 
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Fig. 6. Cloudy predictions of the line intensities as a function of ionization parameter, U, and radiation field. Go- Each connected track correspond 
to a given density (black = 10 1 0 cnT 3 ; light gray = 10 5 0 cirT 3 , in steps of 0.5 dex). The observed ratios are indicated by the yellow (dwarf galaxies) 
and red hashed (normal galaxies) bands. The [O i]/[C n] and [Cii]/Ltir values corrected for non-PDR [Cii] emission. The stars/squares show 
the fiducial Hu region/PDR models for the dwarfs (green symbols) and the normal galaxies (red symbols). Left: results for the high-metallicity 
(2Z ism ) grid of models. Right: results for the low-metallicity (0.25 Z IS m) grid of models. The fiducial model for normal galaxies is added to the 
low-metallicity panel for comparison (smaller red symbols). 


[Nil] 122 line alone. Since we have taken into account multiple 
ionic lines, our derivation of the H n region properties is more 
robust and the previous ionized gas [Cn] fractions may have 
been over-estimated. Alternatively, an unmodeled low-density, 
low-ionization component may increase this fraction by emitting 


more [C ii] (and [N n]), but is not necessary here to match the 
observations. 
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4.3.2. Conditions in the H it region: low-metallicity case 

The observed [Sm] 33 /[Sm]ig ratios are similar in dwarf and 
normal galaxies, and again agree well with densities between 
100 and 1,000cm -3 . From observations, the [Nem]i 5 /[Neii]i 2 , 
[S iv]io/[S hi] is, and [O iii]88 /[N n] 122 ratios are higher in the 
dwarf galaxies than in metal-rich galaxies, which is a clear in¬ 
dication of harder radiation fields (e.g., [Madden et al.|20061 >. At 
low metallicity, the model predictions for those line ratios are 
higher by a factor of ~5 than in the high-metallicity grid (for a 
given set of «h and U) because the radiation field of the star- 
burst is harder. This can be seen by the shift of position of the 
red star from the left to the right panels in Fig. [ 6 ] Even with 
the change of the starburst metallicity, the observations still re¬ 
quire higher U to better match the observed [O m]/[N 11 ] 122 ratio. 
While log U = -3 agrees well with the observed ratios in nor¬ 
mal galaxies, log U = -2.5 agrees better with observations for 
the dwarf galaxies. The [Ne m] is/[Ne 11 ] 12 and [S iv] io/[S 111 ] is ra¬ 
tios are slightly over and under-predicted, respectively, by the 
model, as in the high-metallicity case. Going to higher U, the 
[O iii]/Ltir, which is also higher in the dwarfs than in the normal 
galaxies, is then better reproduced by marginally lower densi¬ 
ties («h — 10 2 cm -3 ), but at the expense of a less well predicted 
[S hi] line ratio. As a best compromise, we take log U = -2.5 
and «h = 10 2 cm -3 as parameters for our fiducial low-metallicity 
H 11 region model. In those conditions, our fiducial model pre¬ 
dicts that only 1 % of the [C 11 ] emission is originating in the ion¬ 
ized gas for the dwarf galaxies. However, the [O iii]/Ltir ratio 
remains under-predicted by all the grid models. We discuss the 
origin of this major discrepancy in Sect. [5] 


4.3.3. Conditions in the PDR: high-metallicity case 


In the PDR, the ratios of [Oi]/[Cii], [Cii]/Ltir, and [Oi]/Ltir 
are sensitive to both density and Gq. With n\\ - 10 3 5 cm -3 and 
logGo — 3.3, the fiducial PDR model for metal-rich galaxies 
matches well all three ratios ( bottom, left panels of Fig. [ 6 }. This 
diagram also shows that, if more [C 11 ] emission were originating 
from the PDR than the 50% assumed here, reducing the density 
and Go values by 0.5 dex would suffice. 

However, the F60/F100 band ratio is over-predicted by a fac¬ 
tor of ~2 for the metal-rich galaxies. Contrary to the dwarfs, a 
substantial fraction of the dust emission at wavelengths >100 gm 
may not be related to the star-forming regions in normal galax¬ 
ies (e.g.. Dale et ak|2007| [Bendo et ah|2012a]|2014||De Looze| 


|et al.|2014b i. If old stars contribute to dust heating, and thus to 

the FIR emission, the line-to-Z/riR ratios associated with the star 
formation would be higher than observed, which would imply a 
lower Go than deduced currently. Alternatively, since F60/F100 
is sensitive to Ay and all models are stopped at Ay of 10 mag, the 
IRAS band ratios of the normal galaxies may be under-estimated 
by the models. 


4.3.4. Conditions in the PDR: low-metallicity case 

Considering the low-metallicity case, the predictions for 
[Oi]/[Cii] decrease by a factor of ~2 ( bottom, right panel of 
Fig. | 6 } owing to a decrease in the PDR temperature ( jRollig et al.| 
|2006p ~ Moreover, in the models with U > 10 -2 , we observe two 
effects. First, those models have low densities and as a conse¬ 
quence the thickness of the cloud becomes larger than the inner 
radius. As a result, there is significant geometric dilution of the 
UV radiation field. The C + and O 0 emitting layers are larger, and 
thus [C ii]/Ltir and [O i]/Ttir ratios are high. Second, regardless 


of metallicity, there is more dust absorption in the H 11 region at 
those high-G values ( Abel et al.|2009j ), and thus little increase of 
Go in the PDR for increasing U. As a consequence, the predicted 
ratios of [Oi]/[Cii], [Cii]/Ltir and [Oi]/Ltir flatten out with U. 

The observed [Oi]/[Cii] ratios in the dwarfs and normal 
galaxies are similar, while [Cii]/Ljir and [Oi]/Ltir are higher 
in the dwarfs. The ratios of the dwarfs indicate lower Go or 
higher n H values. We find that a solution with n H - 10 4 cm -3 
and logGo - 2.5 is adequate to reproduce all three ratios. This 
PDR solution agrees with the observed F60/F100 band ratio for 
the dwarfs. 


5. Discussion 


5.1. Interpretation in terms of ISM structure 

The fiducial models imply different sizes for the low-metallicity 
and high-metallicity components (Table |3j. The ionic lines are 
emitted by a modeled shell, for which the size can be converted 
to an effective Stromgren radius If. eff- which can be thought of 
as a cumulative volume of Hu regions in a galaxy. As an ex¬ 
ample, for the high-metallicity case, Rs e ff of 16 pc corresponds 
to about one hundred H 11 regions if we consider each to have a 
typical radius of 3-4pc. Because the fiducial low-metallicity Hu 
region model has harder radiation field and lower density, we 
find that its effective radius is three times larger than that of the 
high-metallicity H 11 region. The depth of the PDR, defined as 
the location where the gas is half atomic half molecular, is sim¬ 
ilar in the low-metallicity and high-metallicity ISM. However, 
the low-metallicity PDR is located further away. 

Discrepancies in the predicted [O iir|/L- nR ratios between the 
high and low-metallicity grids are mainly due to a reduced gas- 
phase metallicity. This reduces the absolute [0 111 ] emission pro¬ 
duced in the Hii region, while both grids of models have the 
same, fixed input luminosity which is fully reprocessed in the 
PDR, exiting as Ltir (no escaping UV). To increase the predicted 
ratios of ionic lines to L TIR and better match the observations, we 
would need to reduce significantly the mean Ay of the models. 
This can be achieved by accounting for an escape fraction, re¬ 
ducing the covering factor of the PDR (hence L TIR and the aver¬ 
age column density). In other words, the emission from the PDR 
must be reduced relative to that of the H 11 region, by a factor 
of approximately 3 in our case (for the green star to fall on the 
observed median [O iii]/Ltir value in the right panel of Fig. | 6 j. 
Lower PDR covering factors were also required in the detailed, 
multi-phase modeling of Haro 11 (jCormier et a k[2012j >. The fact 
that the median value of Ltir/Eb- which represents the fraction 
of obscured/unobscured emission, is ~2 in the B08 galaxies and 
0.7 in the DGS galaxies corroborates our finding of lower PDR 
covering factors for the low-metallicity models. This also con¬ 
firms the observed trend of decreasing [O m]/Gn R with increas¬ 
ing Lx IR /L b . The high escape fraction in dwarf galaxies is also 
observed in their SEDs (e.g., Galliano et al. |2008[ >. 

The higher ratios of [Oi] and [Cn] to Ltk indicate lower 
U in the low-metallicity PDR, which can be interpreted as a 
change of ISM configuration. Since the Hu region and PDR 
low-metallicity fiducial models have the same Go (see Fig. | 6 j, 
placing the PDR adjacent to the Hii region and increasing the 
density would result in lowering U. In the high-metallicity fidu¬ 
cial model. Go is higher in the PDR than in the H 11 region, which 
would require to place the PDR inside the H 11 region. Combined 
with the fact that the low-metallicity H 11 region fills a larger vol¬ 
ume than the high-metallicity Hii region, this implies that the 
minimum distance of the low-metallicity PDR to the starburst 
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should be 6 times larger than in the high-metallicity setting. The 
increase in the typical PDR distance is a natural consequence of 
a lower covering factor. 

These considerations highlight a change in the ISM structure 
and PDR distribution of low-metallicity galaxies, where the fill¬ 
ing factor of ionized gas relative to the PDRs is higher than in 
metal-rich galaxies. The extent of the Hu region is bigger and 
the UV photons travel larger distance before they reach a PDR 
clump. 


5.2. Interpretation of the [Cu]/L TU< and [Oi]/L T1R trends 

In Sect. |3.4| we observed that the line-to-LxiR ratios decrease 
with increasing Lxir and Lxir/Lb, and that [Cii]/Ltir is anti¬ 
correlated with F60/F100. Using the PDR models, we discuss 
which physical effects dominate in producing these trends. 

Line-to-LxiR ratios have been studied extensively in massive 
galaxies (normal, star-forming galaxies and ULIRGs). They give 
the budget of the gas cooling over the dust cooling. The ([C n] 
+[0 i](, 3 )/T [ ik ratio, a common measure of the photoelectric ef¬ 
ficiency in the PDR, is observed to vary (see Fig. [5]>, with a 


clear deficit of [C ii]/L T ir in ULIRGs ( |Luhman et al.|2003] Diaz 


Santos et al. 2013 i, and FIR lines in general (Gracia-Carpio et al. 
201 1|>, as well as a decrease with FIR color in star-forming galax 


ies (Malhotra et al. 20011 . Several hypotheses to explain the ap¬ 
parent line deficit in galaxies include ( 1 ) more dust screening 
at high ionization parameter (Luhman et al. 2003; Abel et al. 
2009 ); (2) lower photoelectric efficiency due to charged small 
dust grains under intense radiation fields (Malhotra et al. 2001 


|Croxall et al,|2012] >; (3) lower photoelectric efficiency due to 
decrease in the PAH abundance under hard radiation fields (Mad¬ 
den et al. 2006; Rubin et al. 2009); and (4) large dust optical 


depths leading to extinction of the emission lines (Papadopoulos 
et al,|2010|[Rangwala et al.[20TT) . The dwarf galaxies extend the 
relative importance of the gas cooling to high values. Here, we 
aim to investigate which of those interpretations can explain the 
difference from the normal to the dwarf galaxies. 

Because of their low dust content, variations in grain abun¬ 
dances and charging could affect the photoelectric efficiency dra¬ 
matically in the dwarf galaxies. The quantity Lxir sums emis¬ 
sion from several grain populations; from the big grains (pre¬ 
dominantly) to the very small grains and PAHs, which are the 
most important for the photoelectric heating ( |Bakes & Tielens| 
|1994| >. Those grains emit from all ISM phases, and the fraction 
originating in the PDRs themselves is unknown. Star-forming 
dwarf galaxies generally exhibit enhanced small grain abun¬ 
dances ( [Galliano et al.|2003] ), which correlate with H n regions, 
and reduced PAH abundances, which mostly trace the PDR, as 
a result of the hard radiation fields ( [Madden et al.||2006j ). Thus 
those two populations may balance their relative contribution in 
the photoelectric effect. Moreover, the mean free path of UV 
photons in low-metallicity environments is larger, resulting in a 
UV field dilution over larger spatial scales, which reduces grain 
charging (Israel & Maloney |20 11 1 . Sandstrom et al.|( |2012} find 
that, even though their abundance is reduced, PAHs are gener¬ 
ally more neutral in the star-forming regions of the SMC than in 
typical spiral galaxies. Lebouteiller et al. ( 2012] ) examine both 
spatial distributions of ([C n] 157 +[Oi] 63 ) with Lxir and PAH 
emission in the star-forming region NT 1 B of the LMC, and find 
that PAH emission traces best the gas heating as compared to 
Lxir. The net result of those considerations is that, overall, the 
observed ([C 11 ] 157 +[0 i]63)/Lxir ratios are high in the dwarfs 
(even if we subtract a small contribution of the ionized gas to the 
[C 11 ] emission), indicating either relatively efficient photoelec- 
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trie heating or non-photoelectric excitation of these lines (cos¬ 
mic rays. X-rays, etc.). The ratios tend to flatten out to 1% at 
low luminosity, which could indicate a ceiling on the photoelec¬ 
tric efficiency. The influence of more involved changes to the 
modeling (decrease of PAH abundance, grain size distribution, 
or X-ray illumination) will be examined in Paper II. Here we 
conclude that a simple change in geometry, by putting the cloud 
further away, with viable density and illumination conditions can 
suffice to explain the high observed line-to-LxiR ratios. 

[Cii]/Ltir and, to a lesser extent, [Oi]63/Lxir decrease as 
Go increases (Fig. [6]). Two main effects are responsible for that: 
the increased grain charging and increased dust-to-gas opacity 
which cause a larger fraction of the radiation to go into heat¬ 
ing the dust. The dwarfs show high [C ii]/L T ir and [O i]/L T ir be¬ 
cause Go and Go/«h (which traces grain charging) are lower in 
the low-metallicity PDR than found in normal galaxies ( |Mal- 
|hotra et al. 2001). Since [O i]63 increases faster with Go than 
[Cn], [Oi] 63 /Lxir is less affected, except at the highest Go val¬ 
ues. We also note that, [O i]63 being more sensitive to density 
than [Cn], [Oi^/Lxir increases with n H for a given Go- In 
the models, F60/F100 follows immediately from Go- The anti¬ 
correlation between [C ii]/Lxir and the FIR band ratio, which we 
interpreted as an effect of compactness, is evident in Fig. [6] and 
the lack of correlation between [Oi]63/Lxir and FIR color agrees 
with the models. The stronger dependency of [C ii]/Lxir on Go 
compared to [Oi^/Lxir may account for the larger scatter in 
the observed [Cii]/Lxir values. The positive correlation between 
[O i] 63 /[C 11] 157 and FIR color is also reproduced by the models, 
which may be understood as an aging effect of the star-forming 
regions that expand as they age. The models do not show signif¬ 
icant optical depth effects. The fact that the [0 1)145 line would 
become optically thick at higher column densities than [Oi]g 3 , 
and that [0 1)145 and [O 1)^3 behave similarly with Lxir and FIR 
color in the observations also confirms the absence of optical 
depth effects at those physical scales. Nevertheless, given the low 
statistics that we have for the [01 ] 145 line and that [01 ] 145 is un¬ 
correlated with Lxir/Lr, we cannot confidently rule out optical 
depth effects. 

Overall, the low-metallicity PDRs are described by a low- 
Go regime accompanied by a density increase. This low G o/«h 
causes UV field dilution and reduces grain charging ( |Israel &| 
[Maloney] |201 1| ). The proposed changes in ISM structure can 
therefore explain the higher line-to-LxiR ratios found at lower 
Lxir and Lxir/Lb (i.e. increased porosity with lower filling fac¬ 
tor of dense phases), as noted in Sect. |3.4| 


6. Conclusions 

We present observations of the FIR fine-structure cooling lines 
[Cii] 157 pm, [Oi] 63 and 145 pm, [Ohi] 88 pm, [N 111 ] 57 pm, 
and [N 11 ] 122 and 205 pm obtained with the PACS spectrometer 
onboard Herschel in 48 low-metallicity galaxies from the Dwarf 
Galaxy Survey. 

- The brightest FIR line in the dwarfs is [O 111 ], on average 
twice brighter than [Cn], the [Oiii] 8 s/[C 11 ] 157 ratio vary¬ 
ing from 0.5 to 13. The FIR lines are bright relative to 
the TIR luminosity, together accounting for several percents 
of Lxir. [Cii]/Ltir range from 0.04 to 0.7% and ([C 11)157 
+[Oi)63)/Lxir range from 0.2 to 1 . 0 %, which is about twice 
as high as in metal-rich galaxies. The gas cooling is enhanced 
on galaxy-wide scales. When zooming in onto the active re¬ 
gions of the LMC, line-to-LxiR ratios are generally higher 
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than galaxy-integrated values. This matches our interpreta¬ 
tion that the dwarf galaxies are dominated by emission from 
their star-forming regions even on galaxy integrated scales. 

- When compared to the Brauher et al. ( j' 2008 j i ISO sample, 
the dwarf galaxies occupy a completely different parameter 
space than the more metal-rich galaxies in several central 
tracers. In particular, the [Oiii]88/[Nii]i22, [Niii] 57 /[Nii]i 22 , 
[C n] i57/[N ii] 122, and [O iii] 8 s /[0 i]63 line ratios are very high 
in the dwarfs, which we interpret in terms of radiation 
field hardness and phase filling factor, with prominent high- 
excitation ionized regions. In spite of the possible influ¬ 
ence of abundance variations, the [O m]/[N ii] ratio appears 
to be a good tracer of the radiation field hardness. The low 
[N ii] i22/[C ii] 157 ratios indicate that there is no substantial 
contribution of the ionized gas to the [C n] emission (typi¬ 
cally < 15 %). 

- With the help of radiative transfer models, the FIR line ra¬ 
tios are interpreted as a change in the ISM properties. The 
high [Cii]/Ltir and [Om]/LxiR ratios indicate a combina¬ 
tion of moderate FUV fields and low covering factor of the 
PDR. The structure of the low-metallicity ISM is qualita¬ 
tively and quantitatively different than found in metal-rich 
environments, harboring only a low filling factor of dense 
gas surrounding the star formation sites, embedded in a large 
volume filling factor of diffuse gas. In this leaky structure, 
UV photons escape from Hu regions out to large distances. 
Overall, the filling factor of the ISM affected by the star¬ 
forming regions is high. 


Given the prospects for observing FIR lines at high redshift 
(ALMA, NOEMA, CCAT), the [Cn]i57, [O iii]xs (twice brighter 
than [C n] 157), and [ 0 1]63 (as bright as [C11] 157) lines appear as 
promising diagnostics of the low-metallicity ISM conditions. 
Each line account for about 0 . 3 % of the TIR luminosity. [C11] 157 
is a reliable tracer of the PDR, with little contribution from the 
ionized gas. [OmJss traces the high-excitation ionized gas while 
[Oiii]/[Oi] indicates, independently of elemental abundance 
variations, the filling factor of ionized/neutral phases, which we 
find high in the low-metallicity star-forming galaxies. 
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Table 4. Line fluxes of the DGS galaxies. 
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Appendix A: Details of the PACS observing strategy 


Table A.l. Technical details about the DGS observations. 

Notes. ( a ) Coordinates of the map center. ( b ) Observation Identification num¬ 
ber. (c) Observing mode: UN=Unchopped Line Scan, CN=Chopnod Line Scan, 
WS=Wavelength Switching. When left blank, the coordinates, observing mode, 
and map size are identical to the previous line. 


Source Name Coordinates^ OBSID <6) Mode (c) Map size Line 

02000) (") 2 


Haro 11 

0h36m52.50s 

-33d33ml9.0s 

1342199236 

CN 

51x51 

[C II] 157 , [O i] 145, [Om]n 


0h36m52.50s 

-33d33ml9.0s 

1342199237 



[Ol]63> [N IIl]57 , [N II] 122 


0h36m52.36s 

-33d33ml9.7s 

1342234063 


47x47 

[N 11)205 

Haro 2 

10h32m31.88s 

+54d24m03.7s 

1342221392 

CN 

47x47 

[Oi] 63 




1342230081 


71x71 

[Cn] 157 




1342230082 


47x47 

[Oiii] 88 

Haro 3 

10h45m22.41s 

+55d57m37.4s 

1342221890 

CN 

79x79 

[Om] 88 




1342221891 


79x79 

[Oi] 63 




1342221892 


95x95 

[Cn] 157 




1342221893 


47x47 

[N n]i22 




1342221894 


71x71 

[O 1 ] 145 




1342221895 


63x63 

[N iii] 57 

He 2-10 

8h36ml5.18s 

-26d24m33.9s 

1342221971 

CN 

47x47 

[N m] 57 




1342221972 


47x47 

[O i] 145 




1342221973 


53x53 

[OlIl ] 88 




1342221974 


53x53 

[Oi] 63 




1342221975 


51x51 

[Cn] 157 

HS 0017+1055 

0h20m21.40s 

+ Ildl2m21.0s 

1342234065 

CN 

47x47 

[Cn] 157 




1342212588 



[Om] 88 

HS 0052+2536 

0h54m56.00s 

+25d53m23.0s 

1342213134 

CN 

47x47 

[Cn] 157 




1342213135 



[Oi] 63 




1342213136 



[Om] 88 

HS 0822+3542 

8h25m55.52s 

+35d32m32.0s 

1342220753 

CN 

47x47 

[Cn] 157 

HS 1222+3741 

12h24m36.72s 

+37d24m36.6s 

1342232306 

CN 

47x47 

[Cn] 157 




1342199399 



[O iii] 88 

HS 1236+3937 

12h39m20.14s 

+39d21m04.8s 

1342199400 

CN 

47x47 

[Cn] 157 

HS 1304+3529 

13h21ml9.92s 

+32d08m23.0s 

1342199736 

CN 

47x47 

[Cn] 157 


13h06m24.24s 

+35dl3m41.6s 

1342232552 



[Om] 88 

HS 1319+3224 

13h21ml9.71s 

+32d08m25.2s 

1342199227 

CN 

47x47 

[Cn] 157 




1342232551 



[Om] 88 

HS 1330+3651 

13h33m08.28s 

+36d36m33.1s 

1342199734 

CN 

47x47 

[Cn] 157 




1342199735 



[Om] 88 




1342232549 



[Oi] 63 

HS 1442+4250 

14h44ml2.81s 

+42d37m44.0s 

1342208927 

CN 

47x47 

[Cn] 157 

HS 2352+2733 

23h54m56.70s 

+27d49m59.0s 

1342213133 

CN 

47x47 

[Cn] 157 




1342234066 



[Om] 88 

II Zw 40 

5h55m43.10s 

+3d23m20.0s 

1342228249 

CN 

79x79 

[Oi] 63 


5h55m43.10s 

+3d23m20.0s 

1342228250 


71x71 

[Om] 88 


5h55m43.10s 

+3d23m20.0s 

1342228251 


85x85 

[0 i] 145 


5h55m43.00s 

+3d23m31.4s 

1342228252 


47x47 

[N II] 122 


5h55m43.10s 

+3d23m20.0s 

1342228253 


95x95 

[Cn] 157 

IZw 18 

9h34m01.44s 

+55dl4m34.8s 

1342220973 

CN 

47x47 

[Cn] 157 


9h34m02.03s 

+55dl4m28.1s 

1342253757 



[Oi] 63 


9h34m02.03s 

+55dl4m28.1s 

1342253758 



[Om] 88 

IC 10 

0h20m26.10s 

+59dl7ml7.6s 

1342214364 

UN 

81x81 

[Oi] 63 


0h20m52.12s 

+59dl7m42.4s 

1342214365 


47x47 

[N n]i22 


0h20m56.70s 

+59dl7m46.8s 

1342214366 


47x47 

[0 1 ] 145 


0h20m26.10s 

+59dl7ml7.6s 

1342214367 


81x81 

[0 iii] 88 


0h20m48.32s 

+59dl8m09.0s 

1342214368 


115x183 

[Cnlis, 


0h20ml8.70s 

+59dl8m45.2s 

1342214369 


81x47 

[Om] 88 


0h20ml4.00s 

+59dl9m53.0s 

1342214370 


149x81 

[Cnlis, 


0h20ml8.70s 

+59dl8m45.2s 

1342214371 


81x47 

[Oi] 63 


0h20ml3.30s 

+59dl9m57.5s 

1342214372 


115x47 

[Oi] 63 


0h20ml3.30s 

+59dl9m57.5s 

1342214273 


115x47 

[O iii) 88 


0h20ml8.00s 

+59dl8m44.0s 

1342223371 


47x47 

[N n]i22 


0h20m28.25s 

+59dl6m58.0s 

1342223372 


47x47 

[N n]i22 


0h20m28.25s 

+59dl6m58.0s 

1342223373 


47x47 

[0 1 ] 145 


0h20ml8.00s 

+59dl8m44.0s 

1342223374 


47x47 

[0 1 ] 145 
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Table A.l. Technical details about the DGS observations. 


Source Name 

Coordinates® 

(J2000) 

OB SID® 

Mode® 

Map size 
(")2 

Line 

LMC-30Dor 

5h38m35.00s 

-69d05m39.0s 

1342222085 

UN 

143x47 

[Om] 88 , [Cn] 157 


5h38m48.00s 

-69d06m37.0s 

1342222086 


71x47 

[Oi] 63 


5h38m58.00s 

-69d04m43.0s 

1342222087 


143x47 

[Oi] 63 


5h38m58.00s 

-69d04m43.0s 

1342222088 


143x47 

[OlIl] 88 , [C II] 157 


5h38m48.00s 

-69d06m37.0s 

1342222089 


71x47 

[Omfel, [C II] 157 


5h38m56.66s 

-69d04m56.9s 

1342222090 


47x47 

[0 i] 145 


5h38m34.92s 

-69d06m07.0s 

1342222091 


47x47 

[N II] 122 


5h38m45.00s 

-69d05m23.0s 

1342222092 


191x47 

[Oi] 63 


5h38m35.00s 

-69d05m39.0s 

1342222093 


143x47 

[0 i] 145 


5h38m45.00s 

-69d05m23.0s 

1342222094 


191x47 

[OlIl] 88 , [Cll]i57 


5h38m46.10s 

-69d04m58.8s 

1342222095 


47x47 

[N II] 122 


5h38m45.00s 

-69d05m23.0s 

1342222096 


191x47 

[0 l] 145 


5h38m35.00s 

-69d05m39.0s 

1342222097 


143x47 

[Oi] 63 


5h38m38.00s 

-69d06m00.0s 

1342231279 


123x47 

[C II] 157 


5h38m30.00s 

-69d06m07.0s 

1342231280 


123x47 

[Om] 88 


5h38m30.00s 

-69d06m07.0s 

1342231281 


123x47 

[Cn] 157 


5h38m55.00s 

-69d03m49.0s 

1342231282 


85x47 

[Ol] 63 ,[Cll]l57 


5h38m30.00s 

-69d06m07.0s 

1342231283 


123x47 

[Ol] 63 ,[Cll],5 7 


5h38m56.00s 

-69d04m50.0s 

1342231284 


85x47 

[N II] 122 


5h38m40.00s 

-69d04m38.0s 

1342231285 


85x47 

[Om] 88 

LMC-N11A 

4h57ml 6.00s 

-66d23m24.0s 

1342214638 

UN 

47x47 

[C 111,57, [Ol] 63 

LMC-N11B 

4h56m95.14s 

-66d24m23.8s 

1342188940 

ws 

127x127 

[O iii ] 88 


4h56m95.14s 

-66d24m23.8s 

1342188941 


127x127 

[Cn]i 57 , [Ol]6 3 , [0 i] 145 , [N Il]l 22 , [N h]205 , [N IIl]57 


4h56m47.00s 

-66d24m32.0s 

1342219439 

UN 

47x47 

[N n]i22 


4h56m57.20s 

-66d25ml3.0s 

1342225175 


47x47 

[N II] 122, [0 i] 145 

LMC-N11C 

4h57m48.00s 

-66d28m30.0s 

1342221976 

UN 

47x47 

[OlIl] 88 , [C 111,57 


4h57m48.00s 

-66d28m30.0s 

1342221977 


47x47 

[Ol] 63 


4h57m40.30s 

-66d26m59.0s 

1342221978 


81x47 

[Oi] 63 


4h57m50.59s 

-66d29m49.7s 

1342221979 


81x81 

[C II] ,57 


4h57m50.59s 

-66d29m49.7s 

1342221980 


81x81 

[Ol] 63 


4h57m40.30s 

-66d26m59.0s 

1342221981 


81x47 

[OlIl] 88 , [C II] 157 

LMC-N11I 

4h55m44.00s 

-66d34m24.0s 

1342222769 

UN 

115x47 

[Oi] 63 




1342222770 


115x47 

[OlIl] 88 , [C II] 157 

LMC-N158 

5h39ml3.00s 

-69d3 0ml 8.0s 

1342214032 

UN 

81x81 

[OlIl] 88 , [C It] 157 




1342214033 


81x81 

[Ol] 63 

LMC-N159 

5h39m37.53s 

-69d45m27.0s 

1342222075 

UN 

47x47 

[0 i] 145 


5h40m08.01s 

-69d44m53.0s 

1342222076 


47x47 

[0 l] 145 * [N Il]205 


5h39m39.50s 

-69d46ml0.7s 

1342222077 


47x47 

[0 i] ,45 


5h39m38.00s 

-69d45m48.0s 

1342222078 


81x81 

[OlIl] 88 , [C It] 157 


5h39m52.00s 

-69d45ml2.0s 

1342222079 


115x47 

[01]6 3 , [N n]i22 


5h39m38.00s 

-69d45m48.0s 

1342222080 


81x81 

[Oi] 63 


5h39m52.00s 

-69d45ml2.0s 

1342222081 


115x47 

[OlIl] 88 , [C II] 157 


5h40m08.00s 

-69d44m46.0s 

1342222082 


81x81 

[0 i] 6 3 , [N II] 122 


5h40m08.00s 

-69d44m46.0s 

1342222083 


81x81 

[OlIl] 88 , [C II] 157 


5h39m53.70s 

-69d45ml3.0s 

1342222084 


81x81 

[0 l] 145 * [N Il]205 

LMC-N159S 

5h40m00.00s 

-69d50m33.4s 

1342222767 

UN 

81x81 

[Oi] 63 


5h40m00.00s 

-69d50m33.4s 

1342222768 


81x81 

[C II] 157 

LMC-N160 

5h39m38.00s 

-69d39m04.7s 

1342222071 

UN 

47x47 

[0 i] 145 


5h39m41.30s 

-69d38m50.0s 

1342222072 


75x103 

[Oi] 63 


5h39m44.21s 

-69d38m40.2s 

1342222073 


47x47 

[0 i] 145 


5h39m41.30s 

-69d38m50.0s 

1342222074 


76x105 

[Ollllgg, [C It] 157 

Mrk 153 

10h49m05.04s 

+52d20m07.8s 

1342209014 

CN 

47x47 

[Oi] 63 




1342209015 



[C II] 157 




1342232267 

CN 

47x47 

[0 iii ] 88 

Mrk 209 

12h26ml6.02s 

+48d29m36.6s 

1342199423 

CN 

47x47 

[OlIl] 88 , [C It] 157 




1342199424 



[Oi] 63 

Mrk 930 

23h31m58.29s 

+28d56m49.9s 

1342212518 

CN 

47x47 

[Oi] 63 




1342212519 


47x47 

[0 i] 145 




1342212520 


51x51 

[C II] 157 




1342212521 


47x47 

[N II] 122 




1342212522 


53x53 

[Oiii] 88 

Mrk 1089 

5h01m37.76s 

-4dl5m28.4s 

1342217857 

CN 

47x47 

[Oiii] 88 




1342217858 


47x47 

[0 i] ,45 




1342217859 


51x51 

[C II] 157 




1342217860 


51x51 

[N II] 122 


continued. 


Article number, page 22 of 


30 








Cormier et al.: PACS spectroscopy of the Herschel Dwarf Galaxy Survey 


Table A.l. Technical details about the DGS observations. 


Source Name 

Coordinates® 

(J2000) 

OB SID® 

Mode® 

Map size 
(")2 

Line 




1342217861 


53x53 

[Oi] 63 

Mrk 1450 

Ilh38m35.62s 

+57d52m27.2s 

1342222067 

CN 

47x47 

[Oi] 63 




1342222068 


47x47 

[Oiii ] 88 




1342222069 


47x47 

[N n]i22 




1342222070 


51x51 

[C II] 157 

NGC625 

Ih35m06.00s 

-41d26ml0.3s 

1342222216 

CN 

63x63 

[Om] 88 




1342222217 


47x47 

[Oi] 63 




1342222218 


85x85 

[Cn ] 157 

NGC 1140 

2h54m33.20s 

-10d01m50.0s 

1342214034 

CN 

95x95 

[C II] 157 


2h54m33.54s 

-10d01m42.2s 

1342214035 


63x63 

[Oiii ] 88 


2h54m33.54s 

-10d01m42.2s 

1342214036 


47x47 

[ 01] 63 


2h54m33.54s 

-10d01m42.2s 

1342225171 


47x47 

[N n]i22 


2h54m33.54s 

-10d01m42.2s 

1342224588 


47x47 

[0 i] 145 

NGC 1569 

4h30m49.06s 

+64d50m52.6s 

1342225753 

UN 

95x143 

[C II] 157 




1342225754 


95x79 

[ 01] 63 




1342225755 


119x71 

[Om] 88 

NGC 1705 

4h54ml 3.50s 

-53d21m39.8s 

1342222213 

CN 

63x63 

[Om] 88 




1342222214 


47x47 

[Oi] 63 




1342222215 


71x71 

[C II] 157 

NGC 2366 

7h28m42.60s 

+69dl lml9.0s 

1342220604 


47x47 

[Om] 88 




1342220605 


71x71 

[Cn]i 57 




1342220606 


63x63 

[Oi] 63 

NGC 4214 

12hl5m39.17s 

+36dl9m36.8s 

1342187843 

CN 

lllxlll 

[Oi] 63 




1342187844 



[Om] 88 




1342187845 


95x95 

[C II] 157 




1342188034 



[N Il]205 




1342188035 



[0 i] 145 




1342188036 



[N II] 122 

NGC 4449 

12h28ml 1.90s 

+44d05m39.6s 

1342197813 

UN 

119x191 

[C II] 157 




1342197814 



[01] 63 




1342197815 

CN 

47x47 

[C II] 157 




1342197816 



[Ot]a 




1342197817 



[C II] 157 




1342197818 



[01] 63 




1342197819 



[C II] 157 




1342197820 



[Oi] 63 


12h28ml 1.30s 

44d05m30.6s 

1342223138 

UN 

81x47 

[0 i] 145 


12h28ml 1.30s 

44d05m30.6s 

1342223139 

CN 

81x81 

[Om] 88 


12h28ml 1.30s 

44d05m30.6s 

1342223140 

UN 

81x47 

[N II] 122 


12h28ml5.40s 

44d06m57.6s 

1342223141 


81x115 

[01] 63 


12h28ml5.40s 

44d06m57.6s 

1342223142 


81x115 

[Om] 88 


12h28ml5.40s 

44d06m57.6s 

1342223143 


81x115 

[Cn]i 57 

NGC 4861 

12h59m02.34s 

+34d51m34.0s 

1342208901 

CN 

47x47 

[Ol]6 3 , [0 i] 145 , [Nll]i 22 , [Nil] 205' [N IIl]s 7 


12h59m02.34s 

+34d51m34.0s 

1342208902 


85x215 

[C II] 157 


12h59m02.34s 

+34d51m34.0s 

1342208903 


47x47 

[Om] 88 


12h59m00.30s 

+34d50m43.6s 

1342221887 


47x47 

[Ol]6 3 , [0 i] 145 


12h59m00.30s 

+34d50m43.6s 

1342221888 


47x47 

[Om] 88 

NGC 5253 

13h39m55.96s 

-31d38m24.4s 

1342202125 

CN 

78x78 

[Oi] 63 




1342202126 


94x141 

[C II] 157 




1342214026 

CN 

63x63 

[0 Iii ] 88 




1342214027 


47x47 

[N II] 122 




1342214028 


62x47 

[0 i] 145 

NGC 6822 

19h45m05.00s 

-14d43m22.1s 

1342216633 

UN 

47x47 

[C II] 157 


19h44m50.00s 

14d52m46.9s 

1342216634 


71x71 

[01] 63 


19h45m05.00s 

-14d43m22.1s 

1342216635 


47x47 

[Om] 88 


19h44m50.00s 

14d52m46.9s 

1342216636 


77x47 

[C II] 157 


19h44m50.00s 

14d52m46.9s 

1342216637 


77x47 

[Om] 88 


19h44m52.80s 

-14d43m04.9s 

1342216638 


71x71 

[Oiii] 88 , [Cn]i 5 7 


19h44m52.80s 

-14d43m04.9s 

1342216639 


71x71 

[N II] 122 


19h44m52.80s 

-14d43m04.9s 

1342230147 

UN 

81x81 

[Oi] 63 


19h44m52.80s 

-14d43m04.9s 

1342230148 


47x47 

[0 i] 145 


19h45m05.00s 

-14d43m22.1s 

1342230149 


47x47 

[Oi] 63 

Pox 186 

13h25m48.66s 

-Ild36m37.8s 

1342213284 

CN 

47x47 

[Om] 88 




1342234995 



[C II] 157 

SBS 0335-052 

3h37m44.06s 

-5d02m40.2s 

1342214221 

CN 

47x47 

[C II] 157 


continued. 
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Table A.l. Technical details about the DGS observations. 


Source Name 

Coordinates® 

(J2000) 

OB SID® 

Mode® 

Map size 
(")2 

Line 




1342248295 



[Om] s , 




1342249197 



[Oi] 63 

SBS 1159+545 

12h02m02.37s 

+54dl5m49.5s 

1342199228 

CN 

47x47 

[Cniij, 




1342232309 



[Oiii ] 88 




1342232310 



[Oi] 63 

SBS 1211+540 

12hl4m02.48s 

+53d45ml7.4s 

1342199422 

CN 

47x47 

[Cn] 157 

SBS 1249+493 

12h51m52.53s 

+49d03m26.9s 

1342232266 

CN 

47x47 

[Cn] 157 

SBS 1415+437 

14hl7m01.40s 

+43d30m04.5s 

1342199731 

CN 

47x47 

[Oi] 63 




1342199732 



[0 m ] 88 




1342199733 



[Cn] 157 

SBS 1533+574 

15h34ml3.80s 

+57dl7m06.0s 

1342199230 

CN 

47x47 

[Cn] 157 




1342199231 



[Om] 88 




1342199232 



[Oi] 63 




1342199399 



[Om] 88 

SMC-N66 

0h59m08.00s 

-72dl0m46.4s 

1342214029 

UN 

149x47 

[0 IIl ] 88 , [C Il]l5 7 


0h59m08.00s 

-72dl0m46.4s 

1342214030 


149x47 

[Oi] 63 


0h59m09.73s 

-72dl0m03.4s 

1342214031 


47x47 

[Oi] 63 ,[Cn] 157 

Tol 1214-277 

12hl7ml7.09s 

-28d02m32.7s 

1342199408 

CN 

47x47 

[Cn] 157 




1342234059 



[Om] 88 




1342234060 



[Oi] 63 

UGC4483 

8h37m03.00s 

+69d46m54.0s 

1342203684 

CN 

47x47 

[Cn] 157 


8h37m03.00s 

+69d46m54.0s 

1342203685 



[Om] 88 


8h37m03.00s 

+69d46m31.0s 

1342203686 



[Cn]i5 7 , [Ot] 63 

UM 133 

Ih44m41.28s 

+4d53m25.9s 

1342212533 

CN 

47x47 

[Cn] 157 




1342235699 



[Om] 88 




1342235700 



[Oi] 63 

UM311 

Ihl5m32.00s 

-0d51m38.0s 

1342213288 

CN 

71x71 

[Cn]i5 7 




1342213289 


63x63 

[0 iii ] 88 




1342213290 


47x47 

[N n]i22 




1342213291 


63x63 

[ 01] 63 

UM 448 

1 Ih42ml2.40s 

+0d20m02.7s 

1342222199 

CN 

53x53 

[Om] 88 




1342222200 


47x47 

[N n]i22 




1342222201 


51x51 

[Cn] 157 




1342222202 


51x51 

[0 i] 145 




1342222203 


53x53 

[ 01] 63 

UM461 

1 Ih51m33.35s 

-2d22m21.9s 

1342222204 

CN 

47x47 

[Oi] 63 




1342222205 


51x51 

[Cn] 157 




1342222206 


47x47 

[N II] 122 




1342222207 


53x53 

[Om] 88 

VII Zw 403 

1 Ih27m59.90s 

+78d59m39.0s 

1342199286 

CN 

47x47 

[Nil] 122 , [0 l] 145 




1342199287 


63x63 

[Oiii ] 88 




1342199288 


63x63 

[ 01] 63 




1342199289 


71x71 

[C II] 157 
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Appendix B: Spectral maps and line profiles of the Herschel Dwarf Galaxy Survey 

Display example of the PACS data for NGC 1140. The full atlas is available in the online material. For each spectral line, the flux map and spectrum 
of the brightest spatial pixel, which position is indicated by a black cross on the map, are shown. White circles represent the PACS beam sizes. We 
note that the maps are not PSF-deconvolved. 
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Appendix C: Spitzer IRS fluxes of the compact DGS galaxies 

We have extracted IRS data of the compact DGS galaxies, as described in Sect. |2.2.T| Fluxes of the main MIR fine-structure cooling lines are listed 
here (Table [CT] for the low-resolution data and Table |C.2| for the high-resolution data). 

The spectra of the short and long wavelength modules of the IRS were calibrated using a correction based on the spatial extent of the source. 
The correction, which is wavelength-dependent, assumes a spherical symmetry for the source shape to account for the light that falls outside the 
slit. Using this correction, there is no significant flux offsets between the spectra of the short and long wavelength modules for the low-resolution 
data. For the high-resolution data, extended sources are extracted in different apertures which lead to offsets between the spectra of the short and 
long wavelength modules. In those cases, the short wavelength spectra are stitched to the long wavelength spectra. 

We have compared the overall continuum level of the IRS data to the IRAC and MIPS photometry to verify if all emission is recovered by 
our reduction (i.e. total flux for the galaxy). This is generally the case, except for NGC 1569 which is clearly more extended than the IRS slit. 
For reference, we indicate in Table [C~T| the ratio of the observed-to-synthetic MIPS 24yum flux, denoted Rmipsia- The observed photometry is 
from |Bendo et al.[ ( [2012b) > and the synthetic photometry is measured from the IRS data using the MIPS filter profile. This scaling factor should 
be applied to obtain total IRS fluxes. We opt for the 24 pm emission because it arises from active star-forming regions and it is less affected by 
IRS module stitching than the shorter wavelengths. We note that observed-to-synthetic IRAC 8/rm fluxes are further offset by a factor ~ 1.8 in 3 
galaxies with strong PAH emission (He 2-10, Mrk 1089, NGC 1140). This could be due to PAH (and PDR line) emission outside of the IRS slit as 
those galaxies show extended structures or to calibration uncertainties. 

In Table |C2| we also indicate the average ratio between the high-resolution and low-resolution fluxes of the main spectral lines, denoted 
HR/LR. Line fluxes generally agree within 20% around that ratio. Differences are mainly due to the extraction method and source extent. [Nen] 
and the PAH feature at 12.7 pm, as well as [Oiv] and [Fen] at 25.9/rm, are resolved and thus more accurately measured with the high-resolution 
mode. 
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Table C.l. Low-resolution MIR line fluxes in the compact galaxies of the DGS. 
Notes. Fluxes are in KL 18 W nT 2 . Upper limits are given at 2a. 


Source Name Main Observed Spectral Lines Rmipsia 



[Siv] 10.5 pm 

[Nen] 12.8 pm 

[Ne in] 15.6 pm 

[Siii] 18.7 pm 

[Sin] 33.5pm 


Haro 3 

3.44 ± 0.16 (xlO 2 ) 

3.65 ±0.12 (xlO 2 ) 

- 

- 

- 

- 

Haro 11 

5.28 ± 0.12 (xlO 2 ) 

2.86 + 0.09 (xlO 2 ) 

1.26 ±0.14 (xlO 3 ) 

5.54 + 0.68 (xlO 2 ) 

6.56 + 0.81 (xlO 2 ) 

1.13 + 0.05 

He 2-10 

4.25 + 0.77 (xlO 2 ) 

4.17 + 0.24 (xlO 3 ) 

1.50 ±0.16 (xlO 3 ) 

3.24 + 0.29 (xlO 3 ) 

3.79 + 0.52 (xlO 3 ) 

1.13 + 0.05 

HS 0052+2536 

13.50 + 3.07 

8.55+1.32 

14.97 + 4.51 

21.06 + 6.28 

17.27 + 4.13 

1.15 + 0.09 

HS 0822+3542 

5.23 ± 0.47 

<0.50 

4.35 + 1.54 

< 1.84 

< 2.33 

0.88 + 0.10 

HS 1222+3741 

7.80+1.27 

< 1.81 

< 7.89 

< 6.79 

< 26.26 

- 

HS 1304+3529 

11.10+ 1.58 

< 1.63 

16.24 + 4.76 

16.34 + 4.99 

< 25.44 

1.07 + 0.10 

HS 1319+3224 

<6.94 

< 1.09 

< 3.76 

< 7.97 

< 11.17 

- 

HS 1330+3651 

8.49 + 1.10 

<0.85 

16.45 ± 6.49 

< 17.24 

< 16.75 

- 

HS 1442+4250 

6.92 + 1.58 

< 1.12 

- 

- 

- 

- 

II Zw 40 

1.80 + 0.03 (xlO 3 ) 

83.49 + 5.83 

1.48 ±0.04 (xlO 3 ) 

5.79 + 0.44 (xlO 2 ) 

7.81 + 0.63 (xlO 2 ) 

1.18 + 0.05 

IZw 18 

4.34 + 0.38 

<0.50 

7.42+1.40 

4.44+ 1.19 

3.36+1.66 

1.22 + 0.08 

Mrk 153 

18.91 + 2.00 

5.44+1.33 

29.92 ± 2.54 

19.15+1.77 

24.33 ± 8.66 

1.33 + 0.06 

Mrk 209 

79.92 ± 4.46 

< 8.95 

- 

- 

- 

- 

Mrk 930 

1.29 + 0.20 (xlO 2 ) 

30.19 + 4.27 

2.11 ± 0.18 (xlO 2 ) 

76.56+ 8.33 

1.41 ±0.17 (xlO 2 ) 

1.24 + 0.05 

Mrk 1089 

1.48 ±0.10 (xlO 2 ) 

1.80+ 0.08 (xlO 2 ) 

4.71 ±0.29 (xlO 2 ) 

3.05 ±0.18 (xlO 2 ) 

5.27 + 0.83 (xlO 2 ) 

1.27 + 0.06 

Mrk 1450 

69.77 + 3.19 

10.69+ 1.86 

94.68 ± 3.29 

49.60 + 4.48 

62.92 ± 6.03 

1.07 + 0.05 

NGC 1140 

1.81 ±0.10 (xlO 2 ) 

1.59 + 0.07 (xlO 2 ) 

4.78 + 0.20 (xlO 2 ) 

2.47 ±0.19 (xlO 2 ) 

4.36 + 0.96 (xlO 2 ) 

1.26 + 0.05 

NGC 1569 

- 

- 

2.73 ±0.16 (xlO 3 ) 

1.38 + 0.07 (xlO 3 ) 

2.50 + 0.37 (xlO 3 ) 

3.16 + 0.14 

Pox 186 

36.29 ± 3.66 

0.95 ± 0.25 

13.81 + 3.33 

7.37+ 1.25 

< 4.64 

0.86 ± 0.09 

SBS 0335-052 

20.04 + 2.10 

1.35 + 0.56 

17.35 + 2.71 

4.24+ 1.60 

< 5.81 

1.02 + 0.05 

SBS 1159+545 

3.91 ± 1.78 

< 2.51 

< 3.32 

7.06+ 1.96 

< 6.92 

1.10 + 0.13 

SBS 1211+540 

3.10 + 0.66 

< 1.15 

4.94 ± 0.73 

3.44+ 1.32 

< 16.72 

1.25 + 0.41 

SBS 1249+493 

4.74 + 0.91 

<0.57 

9.85 + 3.89 

< 11.72 

< 16.27 

0.85 + 0.17 

SBS 1415+437 

13.95 ± 0.60 

3.10+1.23 

- 

- 

- 

- 

SBS 1533+574 

28.93 + 2.87 

4.28 + 0.51 

49.49 ± 9.09 

14.10+1.53 

25.28 ± 5.94 

- 

Tol 1214-277 

10.53 + 2.41 

< 1.52 

5.89 + 2.48 

< 5.10 

< 7.11 

1.19 + 0.11 

UM 448 

2.21 + 0.11 (xlO 2 ) 

3.76 + 0.22 (xlO 2 ) 

6.93 + 0.83 (xlO 2 ) 

3.24 ±0.19 (xlO 2 ) 

6.63 + 0.48 (xlO 2 ) 

1.03 + 0.04 

UM461 

45.89 ± 4.23 

2.53 ± 0.62 

25.61+6.50 

11.91 ± 1.88 

11.96 + 3.76 

1.03 + 0.10 

Source Name 


Other Observed Spectral Lines 




[Arn] 6.99 /mi 

[Nev] 14.3 pm 

[Oiv] 25.9 pm 

[Sin] 34.8pm [Nem] 36.0pm 


Haro 3 

< 56.13 

- 

- 

- 

- 


Haro 11 

37.83 ± 11.61 

< 91.56 

< 35.97 5.16+ 1.31 (xlO 2 ) < 

3.75 (xlO 2 ) 


He 2-10 

4.94 + 0.43 (xlO 2 ) 

1.72 + 0.73 (xlO 2 ) 

<9.42 (xlO 2 ) 2.41 ±0.42 (XlO 3 ) < 

1.24 (xlO 3 ) 


HS 0052+2536 

< 4.04 

<7.35 

< 3.11 

25.71 ± 7.21 

- 


HS 0822+3542 

< 2.44 

< 2.14 

<0.85 

2.45 ± 0.92 

< 6.18 


HS 1222+3741 

< 3.04 

< 3.99 

< 7.75 

< 53.57 

- 


HS 1304+3529 

< 8.49 

<6.02 

< 6.11 

< 24.22 

< 16.16 


HS 1319+3224 

- 

< 5.44 

< 1.34 

9.70 + 3.12 

< 24.74 


HS 1330+3651 

< 3.38 

< 16.31 

<4.68 

< 25.65 

< 35.09 


HS 1442+4250 

< 4.40 

- 

- 

- 

- 


II Zw 40 

< 34.98 

< 62.44 

96.47 ± 39.45 4.51 ± 0.49 (x 10 2 ) < 

1.99 (xlO 2 ) 


IZw 18 

< 1.36 

<5.05 

<2.96 

8.10 + 2.76 

< 10.18 


Mrk 153 

< 14.23 

< 5.44 

<4.89 

18.98 + 7.45 

< 21.97 


Mrk 209 

< 26.98 

- 

- 

- 

- 


Mrk 930 

< 28.00 

<9.00 

<10.31 1.08 ±0.08 (xlO 2 ) 

< 68.98 


Mrk 1089 

24.62 ± 7.03 

21.44 + 7.14 

< 13.54 2.64 + 0.51 (xlO 2 ) < 

o 

VO 

00 


Mrk 1450 

< 10.97 

< 14.09 

< 8.63 

31.85 + 3.44 

< 11.22 


NGC 1140 

< 21.00 

< 12.66 

< 13.72 3.28 ±0.30 (xlO 2 ) 

< 75.12 


NGC 1569 

- 

29.45 ± 8.33 

<6.17 (xlO 2 ) 1.61 ±0.22 (XlO 3 ) < 

o 

£ 

(T) 

00 


Pox 186 

< 1.51 

<2.36 

<0.36 

< 4.57 

< 7.85 


SBS 0335-052 

< 2.22 

<6.33 

<5.01 

5.61 ± 1.89 

< 11.93 


SBS 1159+545 

< 11.99 

<5.40 

<0.95 

< 7.34 

< 7.62 


SBS 1211+540 

< 3.13 

< 2.75 

< 2.01 

< 10.29 

< 35.40 


SBS 1249+493 

<9.53 

< 20.34 

< 4.54 

< 45.09 

< 27.83 


SBS 1415+437 

< 4.67 

- 

- 

- 

- 


SBS 1533+574 

<2.66 

< 7.76 

<9.59 

20.62 ± 2.84 

< 11.13 


Tol 1214-277 

< 16.92 

<4.59 

<5.61 

<7.28 13.18 + 3.27 


UM 448 

55.15 + 25.71 

< 38.84 

< 10.29 4.89 + 0.49 (xlO 2 ) < 

1.28 (xlO 2 ) 


UM461 

< 4.73 

< 8.43 

< 3.83 

11.87 + 2.84 

< 16.98 
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Table C.2. High-resolution MIR line fluxes in the compact galaxies of the DGS. 

Notes. Fluxes are in 10 _I8 Wnr ; . Upper limits are given at 2a. The symbol * indicates point-like sources reduced with the optimal extraction 
method. The other sources are extracted from the full apertures. We note that the observation of Pox 186 (AOR key 9007360) is not included 
because it is not centered on the source. 


Source Name Main Observed Spectral Lines HR/LR 



[Siv] 10.5/mi 

[Nen] 12.8/rm 

[Nem] 15.6/jm 

[Sin] 18.7 pm 

[Sin] 33.5 pm 


Haro 3 

4.08 ± 0.17 (xlO 2 ) 

3.52 ±0.13 (xlO 2 ) 

9.84 ±0.74 (XlO 2 ) 

5.03 + 0.39 (xlO 2 ) 

8.51 + 0.24 (xlO 2 ) 

1.08 

Haro 11* 

4.94 ±0.11 (xlO 2 ) 

3.27 + 0.09 (xlO 2 ) 

1.12 ±0.05 (XlO 3 ) 

5.31 + 0.29 (xlO 2 ) 

8.17 + 0.70 (xlO 2 ) 

1.04 

He 2-10 

3.27 ± 0.13 (xlO 2 ) 

3.80+ 0.13 (xlO 3 ) 

1.56 + 0.05 (xlO 3 ) 

2.67 ±0.19 (xlO 3 ) 

3.17 ±0.10 (xlO 3 ) 

0.88 

HS 0822+3542* 

6.43 + 0.45 

< 0.24 

2.63 + 0.21 

1.32 + 0.22 

- 

0.92 

HS 1442+4250* 

8.21 + 0.30 

<0.51 

3.83 + 0.68 

1.69 + 0.59 

3.10+1.03 

1.19 

II Zw 40 

2.00 ± 0.10 (xlO 3 ) 

73.52 ±7.93 

1.41 ±0.09 (xlO 3 ) 

5.21 ±0.22 (xlO 2 ) 

7.82 + 0.21 (xlO 2 ) 

0.97 

IZw 18 

5.69 + 1.42 

< 2.91 

7.89+ 1.02 

5.04+1.29 

3.94 + 0.91 

1.17 

Mrk 153 

24.01 ± 5.47 

3.91+0.77 

31.22+ 1.72 

14.19 + 2.32 

28.89 ± 6.97 

0.99 

Mrk 209* 

85.71 + 2.88 

2.45 + 0.31 

54.49 ± 1.74 

24.34 ± 1.49 

49.69 ± 5.70 

1.07 

Mrk 930 

1.39 + 0.05 (xlO 2 ) 

39.39 ± 2.44 

2.12 + 0.06 (XlO 2 ) 

94.43 ± 4.98 

1.28 + 0.11 (xlO 2 ) 

1.12 

Mrk 1450* 

72.84+ 1.22 

9.18 + 0.72 

93.06 ±3.32 

43.88 ± 1.54 

70.17 + 17.05 

0.98 

NGC 1140 

1.80+0.10 (xlO 2 ) 

1.83 + 0.09 (xlO 2 ) 

6.28 + 0.40 (xlO 2 ) 

3.31 ±0.17 (xlO 2 ) 

4.34 + 0.14 (xlO 2 ) 

1.16 

NGC 1569 

2.47 + 0.03 (xlO 3 ) 

3.05 ± 0.27 (xlO 2 ) 

3.24 + 0.13 (xlO 3 ) 

1.31 ±0.05 (xlO 3 ) 

1.85 + 0.06 (xlO 3 ) 

0.96 

SBS 0335-052* 

14.75 + 0.37 

0.72 ± 0.09 

12.37 ± 0.74 

3.49+ 1.05 

4.70+1.58 

0.71 

SBS 1159+545* 

<5.10 

< 2.19 

< 2.97 

< 4.23 

< 12.83 

n.a. 

SBS 1415+437* 

12.02 ± 0.92 

2.31+0.29 

11.33 + 0.99 

7.58+1.20 

14.04 ± 3.42 

0.81 

Tol 1214-277 

8.86 ± 0.47 

< 1.56 

5.57 + 0.55 

< 2.99 

< 11.06 

0.89 

UM 448 

2.32 + 0.04 (xlO 2 ) 

3.14 ±0.12 (xlO 2 ) 

6.38 + 0.25 (xlO 2 ) 

3.21 ±0.12 (xlO 2 ) 

4.92 ±0.21 (xlO 2 ) 

0.92 

UM461* 

46.71 + 0.90 

1.58 + 0.39 

29.11 ± 1.24 

8.71 + 0.80 

12.64 ± 4.60 

0.86 

NGC 625 

5.72 + 0.07 (xlO 2 ) 

1.66 + 0.05 (xlO 2 ) 

7.46 + 0.26 (XlO 2 ) 

4.40 + 0.15 (xlO 2 ) 

7.78 + 0.39 (xlO 2 ) 

n.a. 

UM 311 

36.80 + 2.91 

43.77 ± 13.69 

1.03 + 0.05 (xlO 2 ) 

76.28 ± 5.09 

59.99 ± 14.60 

n.a. 

VII Zw 403 

14.47 ± 1.42 

5.09 ± 0.34 

23.05 ± 2.60 

18.49 + 2.88 

< 13.60 

n.a. 


Source Name Other Observed Spectral Lines 



[Arn] 6.99 pm 

[Nev] 14.3 pm 

[Oiv] 25.9/jm 

[Sin] 34.8yum 

[Nem] 36.0 pm 

Haro 3 

- 

< 14.50 

17.93 + 5.00 

3.99 + 0.24 (xlO 2 ) 

1.17 + 0.29 (xlO 2 ) 

Haro 11* 

- 

< 5.29 

43.90 + 7.84 

5.58 + 0.45 (xlO 2 ) 

< 1.92 (XlO 2 ) 

He 2-10 

- 

< 20.84 

79.15 ±25.14 

1.93 + 0.04 (xlO 3 ) 

< 2.92 (xlO 2 ) 

HS 0822+3542* 

- 

< 1.02 

- 

- 

- 

HS 1442+4250* 

- 

< 0.43 

2.37 ± 0.74 

3.40+ 1.39 

< 7.32 

II Zw 40 

- 

< 7.46 

79.36 ± 11.53 

3.68 + 0.35 (xlO 2 ) 

1.74 + 0.39 (xlO 2 ) 

IZw 18 

- 

< 4.02 

<4.08 

< 4.37 

< 10.87 

Mrk 153 

- 

< 3.45 

5.10+ 1.92 

23.78 + 6.17 

< 15.73 

Mrk 209* 

- 

< 1.87 

13.13 + 2.81 

32.37 ± 4.07 

16.53 ± 7.43 

Mrk 930 

- 

< 3.29 

12.02 + 2.65 

1.27 ±0.12 (xlO 2 ) 

34.56+ 13.98 

Mrk 1450* 

- 

< 1.04 

8.31 ±0.88 

30.72 + 5.42 

< 15.28 

NGC 1140 

- 

< 9.94 

19.84 + 3.91 

2.71+0.15 (xlO 2 ) 

33.81 ± 16.15 

NGC 1569 

- 

< 8.18 

3.19 ±0.15 (xlO 2 ) 

1.15 + 0.03 (xlO 3 ) 

<5.77 (xlO 2 ) 

SBS 0335-052* 

- 

< 1.63 

5.22+ 1.15 

< 10.55 

< 17.39 

SBS 1159+545* 

- 

< 3.64 

< 2.35 

< 13.85 

< 18.92 

SBS 1415+437* 

- 

< 1.17 

6.51 ± 2.02 

17.77 + 3.91 

< 7.60 

Tol 1214-277* 

- 

< 0.75 

5.43 ± 2.08 

< 8.89 

- 

UM 448 

- 

< 5.64 

30.85 ± 6.35 

4.92 + 0.21 (xlO 2 ) 

< 77.59 

UM461* 

- 

< 0.72 

2.77 ± 0.90 

8.76 + 2.71 

< 20.95 

NGC 625 

- 

< 31.40 

< 3.89 

1.85 ±0.12 (xlO 2 ) 

68.04 ± 22.33 

UM 311 

- 

< 2.55 

6.39 ± 2.62 

56.23 + 5.82 

< 17.00 

VII Zw 403 

- 

< 5.13 

< 1.78 

< 12.35 

< 16.00 
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Appendix D: IRAS and PACS photometry of the compact sample 


Table D.l. L T ir, IRAS 60/100, and PACS 70/100 band ratios. 

Notes, n.a not observed in the PACS 70/rm or 100 pm band. L TiR are SED-integrated values between 3-1 lOOyum, in units of solar luminosity 
for the compact sample and W m~ 2 sr~* (peak surface brightness) for the extended sample. The synthetic IRAS fluxes, indicated in brackets, are 
derived from the SED modeling ofjRemy-Ruyer et~aL](in prep/j). 


Source name 

60/100 

70/100 

Etir 

Compact sample 

Haro 11 

1.21 

1.26 

1.98 x 10 u 

Haro 2 

0.88 

0.93 

6.38 x 10 9 

Haro 3 

0.73 

0.92 

5.43 x 10 9 

He 2-10 

0.89 

1.02 

5.49 x 10 9 

HS 0017+1055 

[1.15] 

1.40 

9.70 x 10 s 

HS 0052+2536 

[0.72] 

1.06 

1.62 x 10 10 

HS 0822+3542 

[0.67] 

1.11 

1.05 x 10 7 

HS 1222+3741 

[0.88] 

0.69 

2.40 x 10 9 

HS 1236+3937 

- 

0.83 

- 

HS 1304+3529 

[0.80] 

0.80 

1.69 x 10 9 

HS 1319+3224 

[0.82] 

0.98 

2.37 x 10 s 

HS 1330+3651 

[0.63] 

0.83 

1.37 x 10 9 

HS 1442+4250 

- 

1.68 

- 

HS 2352+2733 

- 

2.52 

- 

II Zw 40 

0.99 

1.11 

2.89 x 10 9 

IZw 18 

[1.74] 

2.46 

3.34 x 10 7 

Mrk 1089 

0.67 

0.94 

3.68 x 10 10 

Mrk 1450 

0.49 

1.22 

3.03 x 10 8 

Mrk 153 

0.59 

0.93 

1.05 x 10 9 

Mrk 209 

0.97 

0.95 

2.98 x 10 7 

Mrk 930 

0.59 

0.81 

1.95 x 10 10 

NGC 1140 

0.68 

0.95 

3.73 x 10 9 

NGC 1569 

0.93 

1.05 

1.27 x 10 9 

NGC 1705 

0.48 

0.80 

6.29 x 10 7 

NGC 2366 

0.75 

0.83 

1.54 x 10 s 

NGC 4214 

0.62 

0.76 

5.34 x 10 s 

NGC 4861 

0.76 

1.02 

3.51 x 10 s 

NGC 5253 

1.00 

1.01 

1.63 x 10 9 

NGC 625 

0.66 

0.61 

2.85 x 10 s 

Pox 186 

[0.73] 

0.74 

5.64 x 10 7 

SBS 0335-052 

[2.76] 

2.30 

1.76 x 10 9 

SBS 1159+545 

[1.03] 

1.01 

2.00 x 10 s 

SBS 1211+540 

[1.28] 

1.91 

2.86 x 10 7 

SBS 1249+493 

[0.75] 

0.93 

9.34 x 10 8 

SBS 1415+437 

[1.00] 

1.15 

6.54 x 10 7 

SBS 1533+574 

0.63 

0.75 

2.04 x 10 9 

Tol 1214-277 

[0.83] 

0.96 

8.84 x 10 8 

UGC4483 

[0.96] 

0.43 

2.60 x 10 6 

UM 133 

[0.73] 

0.97 

1.13 x 10 s 

UM311 

[0.44] 

0.57 

5.61 x 10 9 

UM 448 

0.93 

n.a. 

9.00 x 10 10 

UM461 

[0.93] 

1.48 

7.24 x 10 7 

VII Zw 403 

0.43 

0.84 

1.99 x 10 7 

Extended sample 

IC 10 

[0.93] 

0.68 

7.30 x 10” 3 

LMC-30Dor 

[1.47] 

n.a. 

1.11 x 10- 3 

LMC-N11A 

[0.90] 

n.a. 

7.92 x 10~ 5 

LMC-N11B 

[1.16] 

n.a. 

1.67 x 10~ 4 

LMC-N11C 

[1.09] 

n.a. 

8.99 x 10- 5 

LMC-N111 

[0.90] 

n.a. 

3.01 x 10~ 5 

LMC-N158 

[1.08] 

n.a. 

2.41 x 10~ 4 

LMC-N159 

[1.18] 

n.a. 

4.79 x 10~ 4 

LMC-N160 

[1.34] 

n.a. 

9.34 x 10~ 4 

NGC 4449 

[1.06] 

0.65 

2.65 x 10~ 5 

NGC 6822 

[0.87] 

0.86 

2.78 x 10~ 5 

SMC-N66 

[1.15] 

n.a. 

4.13 x 10~ 5 
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Appendix E: Line broadening and rotation in the PACS maps 

We find that the [O iii] 88 and [O i]® lines, which are bright with smallest instrumental FWHM, are the most often broadened. The line widths are 
typically ~ 100 km s' 1 larger than the instrumental profile in the brightest sources. For the fainter lines, it is more difficult to determine confidently 
whether or not broadening is present. Broadening of the [C n] 157 line can only be accurately calculated for intrinsic line widths larger than 
150 km s _1 because of its large instrumental FWHM of 240 km s -1 . When observed, [N m ] 5 7 is similarly broadened as are [O iii ] 88 and [Oi] 63 . Such 
broadening indicates the possible presence of strong turbulent motions, or that multiple components contribute to the emission of those lines. The 
fact that we find similar intrinsic line widths for all tracers, despite their different origins, indicates that either they are spatially mixed within our 
beam (which is certainly the case for our compact sources), or that they are affected similarly by the source of this broadening (internal kinematics). 
Several galaxies also display clear signs of rotation in all observed spectral lines. 


Table E.l. Broadening of the FIR lines (in kins 1 ). 


Broadening Rotation 



[N hi] 57 

[Oi] 63 

[O III] 88 

[N 11 ] 122 

[Oi] 145 

[C II] !57 

Axis 

^rotation 

PACS FWHM 

105 

85 

125 

290 

260 

240 



Haro 2 

n. a. 

90 

140 

n. a. 

n. a. 

- 

- 

- 

Haro 3 

100 

90 

110 

- 

- 

- 

NEE-SWW 

±30kms~' 

Haro 11 

160 

130 

200 

200 

100 

160 

N-S 

±50kms~ 1 

He 2-10 

160 

100 

170 

- 

- 

- 

NE-SW 

±20kms~ 1 

HS 1330+3651 

n.a. 

- 

150 

n.a. 

n.a. 

- 

- 

- 

II Zw 40 

n. a. 

80 

110 

- 

- 

- 

N-S 

±30kms~‘ 

Mrk 1089 

n. a. 

80 

130 

- 

- 

- 

NNE-SSW 

±80kms~' 

Mrk 153 

n. a. 

80 

120 

n. a. 

n. a. 

- 

NE-SW 

±30kms~ 1 

Mrk 930 

n. a. 

120 

120 

- 

- 

- 

N-S 

±30kms~‘ 

NGC 1140 

n. a. 

100 

100 

- 

- 

- 

NE-SW 

±40kms~ 1 

NGC 1569 

n. a. 

80 

100 

n. a. 

n. a. 

- 

- 

- 

NGC 1705 

n. a. 

- 

100 

n. a. 

n. a. 

- 

- 

- 

NGC 2366 

n. a. 

80 

70 

n. a. 

n. a. 

- 

- 

- 

NGC 4214 

n. a. 

80 

80 

- 

- 

- 

- 

- 

SBS 1533+574 

n. a. 

- 

80 

n. a. 

n. a. 

- 

- 

- 

UM 448 

n. a. 

150 

150 

- 

80 

150 

NNE-SSW 

±60kms~ l 

IC10 

n. a. 

90 

70 

- 

- 

- 

- 

- 

LMC 

150 

150 

150 

100 

100 

- 

- 

- 

SMC 

n. a. 

150 

150 

n. a. 

n. a. 

- 

- 

- 

NGC 4449 

n. a. 

90 

90 

- 

- 

- 

- 

- 

NGC 5253 

n. a. 

60 

90 

- 

- 

- 

NE-SW 

±30kms~ 1 

NGC 6822 

n. a. 

100 

80 

- 

- 

- 

- 

- 


The broadening (or intrinsic line width) is given as the observed FWHM subtracted by the instrumental FWHM. Uncertainties are on the order of 
20%. We indicate the broadening of the spaxel where the emission peaks. This value can be higher elsewhere because of rotation or instrumental 
effects (uncentered point-source). 
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